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1.  INTRODUCTION 

The  following  paragraphs  provide  background  information  and  objectives 
for  the  ASCOT  effort  as  delineated  in  the  Statement  of  Vtork  (SOW).  Also 
included  is  a  brief  description  of  the  contents  of  the  body  of  the  report. 

1.1  BACKGROUND 

Several  spaceborne  optical  surveillance  and  weapon  system  concepts  of 
interest  to  USAF  and  DARPA  require  precision  Line-of-Sight  (LOS)  and  figure 
control  in  order  to  achieve  their  missions.  To  date,  three  major  system  con¬ 
cepts  have  been  identified:  High  Altitude  Large  Optics  (HALO),  Multi  Mission 
Beam  Control  System  (MMBCS),  and  Large  Optics  Demon  stration  (LODE).  Each 
of  these  system  has  common  features  in  terms  of  large  size,  extreme  com- 
I  plexity,  low  mass,  low  stiffness  and  precision  structural  tolerances.  Each 

concept  is  also  required  to  maintain  stringent  LOS  and  figure  stability 
despite  strong  environmental  and  onboard  disturbances.  The  vibration  dis¬ 
turbance  isolation,  damping,  and  control  technology  required  to  achieve  the 
desired  performance  is  beyond  the  state-of-the-art  practiced  today. 

This  motivates  the  development  and  demonstration  of  new  control  tech¬ 
niques  which  have  broad  application  to  future  large  flexible  military  space 
structures.  Control  system  design  techniques  to  date  have  assumed  comprehen¬ 
sive  and  complex  computer  implementations.  Recognizing  the  ultimate  imple¬ 
mentation  of  these  techniques  in  onboard  space-qualified  digital  computers 
provides  important  practical  reasons  for  their  simplification. 

The  development  and  demonstration  of  such  broadly  applicable,  yet  sim¬ 
plified  and  practical,  techniques  is  the  goal  of  Advanced  Structural  Control 
Techniques  (ASCOT). 
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1.2  OBJECTIVES 


The  objective  of  this  effort  is,  in  general,  the  development,  evaluation, 
and  preparation  for  testing  and  implementation  of  Advanced  Structural  Control 
Techniques  (ASCOT)  that  are  applicable  to  future  large  flexible  military 
spacecraft.  The  two  particular  objectives  are:  (1)  the  development  and 
integration  of  a  Simplified  Systematic  Digital  Design  (S2D2)  which  meets 
anticipated  structural  and  attitude  control  system  performance  requirements, 
and  (2)  the  development  of  an  ASCOT  Technical /Management  Program  Development 
Plan  and  selection,  evaluation,  and  comparison  of  candidate  control  system 
design  techniques  including  S202. 

1.3  DESCRIPTION  OF  REPORT  BODY 

Pursuant  to  the  first  objective  stated  above.  Task  1  includes  those 
activities  related  to  digital  design  technique  development  and  is  reported  on 
in  Section  2.1.  Section  2.2  concerns  Task  2  -  test  planning  and  coordination, 
test  support  and  control  technique  selection  -  and  answers  the  objective  set 
forth  in  Section  1.2.2.  Because  the  tasks  and  subtasks  were  conceived  based 
upon  these  objectives  rather  than  the  input/output  relationship  between 
tasks,  at  certain  points  in  the  text  the  reader  will  be  referred  to  later 
report  sections  for  additional  information. 

Section  2.1.1  specifies  a  set  of  intrinsic  performance  goals  for  the 
S2D2  controller  and  design  technique  in  the  form  of  a  list  of  desired  quali¬ 
ties.  Section  2.1.2  then  specifies  three  typical  disturbances  to  be  applied 
in  evaluating  the  performance  of  S2D2. 

Synthesis  of  the  candidate  S2D2  controller  is  delineated  in  Section 

2.1.3  and  is  followed  by  application  of  the  candidate  technique  to  the 
ASCOT-I  structural  model  described  in  Section  2.2.2. 


Attention  Is  then  turned  to  the  design  of  prefi Iters  (both  analog  and 
digital)  In  Section  2.1.5.  Aliasing  and  the  use  of  prefi Iters  are  discussed 
In  general,  and  the  non-uniform  sampling  period  digital  prefllter  Is  pursued 
particularly.  Both  theoretical  background  and  simulation  results  for  the 
digital  and  analog  prefllters  are  provided. 

Section  2.1.6  provides  post  analysis  studies  of  the  system  designed 
under  Section  2.1.4  In  both  the  frequency  and  time  domains.  Performance  Is 
evaluated  with  respect  to  the  evaluation  criteria  discussed  In  Section  2.2.3. 

A  study  of  which  spacecraft  parameters  are  likely  to  affect  the  control¬ 
ler  performance  Is  given  in  2.1.7  with  the  result  being  a  "perturbed"  version 
of  the  "RTOP"  model  procured  under  2.2.5  which  Is  used  In  the  robustness 
studies  of  Section  2.1.8.  These  studies  are  conducted  using  the  digital  con¬ 
troller  described  In  Section  2.2.5  along  with  the  perturbed  model. 

Description  of  the  Task  2  effort  gets  underway  with  a  brief  account  of 
the  development  of  the  Technical /Management  Program  Development  Plan  which  was 
delivered  midway  of  the  contract  period  as  a  stand-alone  document. 

The  content  of  Sections  2.2.2  and  2.2.3  has  been  discussed  In  the  pre¬ 
ceding  paragraphs.  Section  2.2.4  Includes  a  discussion  of  three  LSS  control 
system  design  technique  candidates  other  than  S2D2. 

Finally,  Section  2.2.5  covers  the  application  of  the  S2D2  candidate  to 
the  RTOP  model  which  Is  a  more  advanced  structure  than  the  ASCOT-I  model. 


2.  TECHNICAL  ACTIVITIES 


2.1  TASK  1  -  DIGITAL  DESIGN  TECHNIQUE  DEVELOPMENT 

2.1.1  Specification  of  Performance  Goals 

The  following  areas  of  importance  have  been  identified  for  practical 
design  techniques  applicable  to  large  flexible  military  spacecraft.  These 
are  the  goals  of  the  S2D2  technique. 

•  Simplicity  of  Controller 

A  digital  controller  design  technique  for  LSS  must  lead  to  controllers 
of  reasonable  order  despite  the  large  order  system  models  which  are  likely  to 
be  used  in  the  design  of  controllers  for  flexible  spacecraft.  This  is  to 
minimize  the  computational  burden  of  the  onboard  computer  when  the  controller 
is  implemented.  This  also  facilitates  simpler  checkout  procedures  for  flight 
qualification  of  control  software. 

•  Straightforward  and  Traceable  Design  Procedure 

The  S2D2  design  technique  should  be  readily  understandable  and  the 
design  process  should  trace  effects  of  closed  loop  control  upon  the  system 
behavior  throughout  the  design  process.  In  this  way  a  designer  can  see  how 
the  system  is  evolving  during  the  design  process,  and  therefore  have  insight 
into  problems  and/or  causes  of  problems  should  they  arise. 

•  Stability  of  Closed  Loop  System 

The  S2D2  technique  should  inherently  provide  stability  of  the  closed 
loop  system  resulting  from  the  combination  of  the  digital  controller  and  the 
LSS  model.  In  fact,  a  reasonable  amount  of  relative  stability  should  be 
inherent. 


•  Inherent  Robustness  Checks 

Robustness  ts  of  particular  concern  to  the  LSS  control  system  designer 
because  accurate  models  are  not  anticipated.  Hence,  the  S202  technique 
should  produce  designs  with  reasonable  robustness  with  respect  to  model  inac¬ 
curacies  and  plant  variation.  In  addition,  the  design  technique  should  pro¬ 
vide  built-in  checks  for  robustness  at  stages  in  the  design  process.  This  is 
also  a  part  of  an  easily  traceable  design  procedure. 

•  Disturbance  Rejection 

Disturbance  rejection  is  a  major  concern  in  control  system  design  for 
LSS  and  should  inherently  be  achieved  through  the  design  process. 

•  Digital  Design  Accomplished  in  Digital  Domain 

Design  of  digital  control  systems  should  be  accomplished  in  the  digital 
domain  so  that  the  effects  of  sampling  and  computational  transport  lags  can 
be  accounted  for  during  the  design  phase.  This  is  in  lieu  of  designing  a 
continuous  controller  and  then  attempting  to  implement  a  digital  equivalent 
which,  at  best,  is  an  approximation  to  the  desired  controller. 

•  Efficiency  of  Design 

The  design  technique  should  be  reasonably  efficient  with  respect  to  com¬ 
puter  processing  and  storage  requirements  and  should  be  algorithmic  in  nature 
so  that  the  design  process  can  be  easily  repeated  as  model  updates  are 
obtained. 

2.1.2  Specification  of  Disturbances 

In  order  to  evaluate  the  noise  rejection  performance  of  the  S2D2  candi¬ 
date  as  applied  to  the  ASCOT-I  and  RTOP  structural  models,  suitable  distur¬ 
bance  signals  must  be  identified.  Of  the  three  distrubance  input  signals 
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applied  In  these  studies,  one  is  a  broadband  random  disturbance  and  two  are 
deterministic  in  nature. 

The  broadband  random  disturbance  is  represented  by  its  Power  Spectral 
Density  (PSD)  as  shown  in  Fig.  1,  and  is  similar  to  that  developed  by 
Riverside  Research  Institute  under  the  ACOSS  program.  The  random  disturbance 
represents  equipment  noise  that  would  be  expected  on  a  "free-flyer"  satel¬ 
lite,  i.e.,  vibrations  of  rotating  machinery  and  fluid  flow.  It  is,  of 
course,  scaled  down  in  magnitude  somewhat  from  that  of  the  ACOSS  application 
because  of  the  considerably  less  massive  spacecraft  model  to  which  it  is 
applied. 

The  two  deterministic  disturbance  inputs  shown  in  Figs.  2  and  3  are 
representative  of  shuttle  orbiter  type  disturbances.  Fig.  2  represents 
astronaut  push-off  within  the  shuttle  orbiter  and  Fig.  3  represents  Reaction 
Control  System  (RCS)  thruster  firing. 

2.1.3  Synthesis  of  Candidate  Primary  Digital  Controller 


The  candidate  S2D2  technique  presented  herein  is  that  of  "ICAT,"  (One 
Controller  At  a  Time).  The  following  sections  delineate  the  theoretical 
background  of  ICAT  and  a  particular  numerical  implementation  of  ICAT. 

Theoretical  Background  of  ICAT.  ICAT  is  a  mul ti-input/mul tioutput  (MIMO) 


digital  controller  design  technique  which  finds  its  basis  in  the  fundamental 
principles  of  classical  analysis  and  design  control  theory.  It  springs  from 
the  fact  that  the  combination  of  a  MIMO  system,  such  as  a  LSS,  and  its  con¬ 
troller  can  be  viewed  as  a  coupled  multiloop  system.  The  controllers  for  the 
loops  cannot  be  designed  independently,  but  they  can  be  designed  one  at  a 
time.  This  is  the  thesis  of  ICAT:  "One  Controller  At  a  Time." 


To  delineate  the  process,  consider  a  system  having  three  inputs  and 
three  outputs.  With  all  possible  feedback  paths  open,  the  transfer  charac¬ 
teristic  between  a  particular  input-output  pair  may  be  examined  and  a 
controller  designed  to  stabilize  the  loop  to  satisfactory  specifications. 
With  this  loop  closed,  another  input-output  pair  may  be  examined  and  a  con¬ 
troller  likewise  designed.  The  second  controller  is  not  designed  indepen¬ 
dently,  because  the  effects  of  the  first  controller  are  taken  into  account 
when  the  first  loop  is  closed.  A  third  controller  may  then  be  designed  with 
the  first  two  loops  closed,  and  so  on  until  all  desired  feedback  paths  are 
closed  through  their  appropriate  controllers.  A  designer  would  probably  not 
choose  to  close  all  nine  possible  feedback  paths  of  this  system;  however,  the 
ICAT  technique  does  not  preclude  this  possibility. 

Three  pertinent  facts  regarding  the  ICAT  technique  and  its  applicability 
bear  mention  at  this  point. 

•  If  the  plant  is  stabil izable,  the  resulting  closed  loop  system  will 
be  stable. 

Stabil izabil ity  simply  implies  that  if  there  are  uncontrollable  and/or 
unobservable  modes,  i.e.,  modes  whose  eigenvalues  cannot  be  changed  by  feed¬ 
back,  their  eigenvalues  must  have  negative  real  parts.  In  this  case  the 
uncontrollable  modes  cannot  result  in  instability,  but  the  controllable  modes 
can.  However,  the  ICAT  approach  can  be  applied  so  that  no  controllable  mode 
can  cause  a  stability  problem  and,  in  fact,  can  produce  a  design  that  will 
ensure  a  specified  amount  of  relative  stability. 

For  example,  suppose  that  the  ICAT  approach  is  applied  to  a  three-loop 
example.  With  all  loops  open,  assume  that  a  controller  is  designed  for  the 
first  loop  so  that  all  controllable  modes  have  closed  loop  eigenvalues  with 
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real  parts  less  than  -a.  Now  with  the  first  loop  closed,  a  controller  Is 
designed  for  the  second  loop  so  that  all  controllable  modes  have  closed  loop 
eigenvalues  with  real  parts  less  than  -a.  Then  with  the  first  two  loops 
closed,  a  controller  is  designed  for  the  third  loop  so  that  all  controllable 
modes  have  closed  loop  eigenvalues  with  real  parts  less  than  -a.  Now  suppose 
a  root  locus  study  is  performed  on  the  first  loop  with  the  second  and  third 
loops  closed  and  with  the  controllers  designed  for  each  loop  included.  The 
controllable  modes  in  this  loop  can  be  separated  into  two  categories:  (1) 
those  modes  controllable  only  from  loop  one  and  (2)  those  modes  controllable 
from  loops  two  and/or  three.  Using  the  gain  factor  for  which  loop  one  was 
designed,  those  modes  controllable  only  from  loop  one  must  have  eigenvalues 
with  real  parts  less  than  -a.,  since  these  modes  are  not  affected  by  the 
designs  in  loop  two  and/or  three,  and  loop  one  was  designed  to  achieve  this 
specification.  The  other  controllable  modes  in  loop  one  must  have  eigen¬ 
values  with  real  parts  less  than  -a  since  these  eigenvalues  are  controllable 
from  loop  two  and/or  three,  which  were  also  designed  to  meet  this  specifica¬ 
tion. 

The  bottom  line  is  that,  as  subsequent  loop  closures  are  made,  eigen¬ 
values  of  preceding  loops  cannot  have  real  parts  greater  than  -a.  However, 
note  that  if  a  subsequent  loop  is  designed  with  a  more  relaxed  specification, 
the  relative  stability  of  the  preceding  loops  can  relax  as  well. 

Although  the  above  arguments  have  been  made  for  the  first  loop  of  a 
three-loop  example,  they  obviously  can  be  extended  to  the  design  of  a  system 
with  many  loops  and  to  loops  other  than  the  first.  In  addition,  other  mea¬ 
surements  of  relative  stability  can  be  used,  e.g.,  gain  margins  and  phase 
margins.  In  order  to  ensure  no  degradation  in  the  relative  stability  of 


loops  previously  closed,  subsequent  loops  should  be  designed  so  that  relative 
stability  is  improved  or,  as  a  minimum,  not  allowed  to  degrade. 

•  Each  controller  in  each  loop  can  be  designed  by  classical  techniques. 

Because  of  the  "one  at  a  time"  nature  of  ICAT,  classical  design  tech¬ 
niques,  such  as  frequency  response,  root  locus,  or  parameter  space  techniques, 
can  be  used  to  design  individual  control  loops.  This  also  makes  possible  the 
use  of  many  tried  and  tested  stability  criteria,  e.g.,  gain  margins,  phase 
margins,  maximum  damping  of  any  mode,  etc.  In  addition,  analysis  and  design 
tools  already  well  developed  are  not  made  obsolete  by  ICAT  but,  to  the  con¬ 
trary,  are  given  a  new  arena  in  which  to  work. 

•  The  effects  of  closure  of  other  loops  upon  a  particular  loop  are 
obvious  by  comparing  the  design  criteria  before  and  after  the 
closures. 

This  intermediate  analysis  of  the  partially  designed  system  gives  the 
designer  insight  into  the  controller  design  problem  at  steps  throughout  the 
design  process.  Particular  problem  areas  may  be  detected  before  the  design 
is  complete  and,  therefore,  attacked  at  their  source.  More  importantly,  the 
designer  is  in  a  better  position  to  discover  what  characteristics  of  his 
design  (or  the  plant  to  be  controlled)  are  causing  stability  and/or  perfor¬ 
mance  problems  than  if  he  uses  a  one-step  design  algorithm. 

Numerical  Implementation  of  ICAT  for  Digital  Controller  Design  Using 
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Frequency  Domain  Techniques.  As  previously  mentioned,  various  design  tech¬ 
niques  are  available  for  accomplishing  a  ICAT  design  of  MIMO  systems.  In 
this  section  the  numerical  implementation  of  ICAT  for  performing  digital  con¬ 
troller  designs  using  frequency  domain  techniques  is  delineated.  The  choice 


of  digital  control  is  made  based  on  the  needed  flexibility  and  expected  hard¬ 
ware  on  systems  for  which  S2D2  is  targeted. 


The  choice  of  frequency  response  techniques  is  based  primarily  on 
obtaining  numerical  efficiency  and  avoiding  numerical  problems.  In  particular 
the  MIMO  systems  for  which  the  ICAT  approach  is  applicable  are  very  high 
order;  thus,  if  root  locus  or  parameter  space  techniques  are  selected,  it 
will  be  necessary  to  determine  z- transforms  of  high  order  systems  and,  con¬ 
sequently,  be  posed  with  numerically  determining  poles  and  zeros  of  z-trans- 
fer  functions.  With  sample  rates  as  expected  for  such  applications,  most  of 
these  poles  and  zeros  will  lie  in  close  proximity  to  z  =  1  in  the  z-plane. 
This  poses  a  significant  problem  in  the  numerical  determination  of  thes^ 
poles  and  zeros,  since  their  locations  will  be  very  sensitive  to  coefficient 
roundoff  in  their  respective  z- transfer  functions.  However,  by  selection  of 
a  frequency  domain  approach  the  tedious  and  numerically  troublesome  z- trans¬ 
form  process  can  be  completely  avoided  as  shown  later. 

Another  need  in  the  numerical  implementation  of  the  ICAT  approach  for 
designing  MIMO  systems  is  a  straightforward  technique  for  the  opening  and 
closing  of  loops  and  the  inclusion  of  controllers  as  the  design  process  pro¬ 
ceeds.  This  is  accomplished  by  modeling  the  MIMO  system  with  vector  -  matrix 
block  diagrams.  In  particular  signals  are  concatenated  to  form  vectors  of 
signals  and  transfer  functions  are  concatenated  to  form  transfer  function 
matrices. 
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The  vector  matrix  model  eases  the  notational  problems  associated  with 
the  handling  of  MIMO  systems  in  the  frequency  domain.  A  transfer  function 
matrix,  [T{s)],  is  defined  to  have  elements,  T^jis),  such  that  Tij(s)  is 


the  transfer  function  or  frequency  response  from  the  jth  input  to  the  ith 
output  of  the  MIMO  system  described  by 

C(s)  =  [T(s)]  R(s)  (2. 1.3-1) 

where  R(s)  and  C(s)  are,  respectively,  the  Laplace  Transforms  or  frequency 
spectrums  of  the  input  and  output  signal  vectors. 

Using  this  notation,  a  block  diagram  of  a  MIMO  system  utilizing  feedback 

and  digital  control  is  shown  in  Fig.  4  where  the  following  definitions  apply: 

n  number  of  modes  used  in  the  design 

analysis 

m  number  of  inputs  and  outputs  of  the 

MIMO  system  under  study 
R(2)  m  X  1  digital  input  vector 

C(s)  m  X  1  continuous  output  vector 

C*(s)  m  X  1  sampled-data  output  vector 

[Di(z)]  m  X  m  digital  feedback  controller  transfer 
function  matrix 

[D2(z)]  m  X  m  digital  forward  path  controller 
transfer  function  matrix 
Im  m  X  m  identity  matrix 

Gho(s)  scalar  zero  order  hold  device  (1  -  e"ST)/s 
T  scalar  sampling  period 

[K]  n  X  m  actuator  modal  gain  matrix 

CGni(s)]  n  X  n  diagonal  modal  transfer  function  matrix 

[L]  n  X  m  sensor  modal  gain  matrix. 

This  block  diagram  has  been  customized  to  the  LSS  design  problem  by  the 
use  of  the  matrices,  [L],  [K],  and  CGni(s)]  to  describe  the  plant,  rather 
than  a  single  m  x  m  matrix,  [T(s)].  Use  of  the  factored  form  allows  the 
modal  gains  to  be  changed  easily  and  efficiently  without  recomputation  of  the 
modal  frequency  responses  in  [G^Cs)]  (exercise  of  (2. 1.3-3)  below).  Another 
assumption  of  this  implementation  is  that  of  an  equal  number  of  inputs  and 
outputs.  This  is  not  required  by  the  ICAT  technique,  and  neither  does  it 
restrict  the  implementation,  because  "dummy"  inputs  and  outputs  may  be 
specified  as  required  at  little  cost  with  regard  to  efficiency.  In  addition, 


this  implementation  is  geared  toward  digital  control;  however,  a  continuous 
control  version  may  be  generated  in  a  similar  manner. 

The  first  step  in  the  ICAT  design  process  is  the  casting  of  the  control 
system  into  the  form  represented  by  the  diagram  of  Fig.  4.  Virtually  any 
MIMO  control  system  can  readily  be  cast  into  this  form  by  determining  what 
the  system  inputs  and  outputs  will  be  and  defining  the  transfer  characteris¬ 
tics  between  them.  For  the  LSS  problem,  using  a  modal  analysis  model  and 
ideal  sensors  and  actuators,  [6n,(s)]  is  diagonal,  and  the  modal  actuator  and 
sensor  gains  are  placed  in  [K]  and  [L],  respectively.  In  the  case  of  sensor 
and  actuator  dynamics,  the  actuator  and  sensor  transfer  functions  may  be 
placed  in  [K]  and  [L]  so  that  no  generality  is  lost.  It  is  easily  seen  that 
in  the  case  where  [K]  and  [L]  are  mth  order  identity  matrices  and  [Gm(s)]  is 
m  X  m  and  not  necessarily  diagonal,  this  Implementation  becomes  the  general 
one  where  [Gm(s)]  =  [T(s)]. 

The  sampled-data  output  vector  of  the  vector-matrix  block  diagram  of 


Fig.  4  is 


C*(s)  =  [In,  +  (LTGn,{s)KGho(s))*D2(z)Di(z)]-l 
[(LTGf,{s)KGho(s))*D2{z)R{z)] 


(2. 1.3-2) 


where  the  symbol,  indicates  a  sampled-data  quantity,  often  called  a 

starred  quantity.  It  should  be  noted  that  a  of  a  matrix  is  a  matrix  in 
which  each  element  is  "*'"ed.  Equation  (1-2)  is  computed  at  discrete  fre¬ 
quency  points  over  the  range  of  interest  of  the  analysis,  with  each  of  the 
transfer  functions  represented  by  its  frequency  response  data  {real  and 
imaginary  parts  vs.  frequency). 

In  evaluating  the  frequency  responses  of  Equation  (2. 1.3-2),  two  types 
of  frequency  responses  must  be  generated  -  those  of  starred  transfer  functions 
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and  those  of  z-transfer  functions.  The  starred  frequency  responses  can 
be  obtained  by  taking  the  z- transform  of  the  quantities  operated  upon  by  the 
starring  operation  and  evaluating  at  z  =  ej“T  or  more  easily  from  the  follow¬ 
ing  Infinite  series: 

1 

G*  (joi)  =  -  I  G(j(d  +  jnus)  +  g{o‘^)/2 

T  n=—  (2. 1.3-3) 

where 

G(s)  Is  the  quantity  being  starred 

T  Is  the  sample  period  In  seconds 

Us  Is  the  sample  rate  In  rad/s 

and 

g(o''‘)  Is  the  Impulse  response  of  G(s)  evaluated  at  t  »  o*^,  I.e., 

11m  [s  G(s)]. 

s+- 

Provlded  that  G(s)  has  more  poles  than  zeros,  the  series  In  Equation 
(2.1.3>2)  converges.  Pragmatic  use  of  (2. 1.3-3)  Is  accomplished  by  trunca¬ 
ting  the  series  so  that  the  frequency  response  is  accurate  to  within  a  spe¬ 
cified  tolerance.  Past  experience  has  shown  that  no  more  than  the  first  20 
positive  and  negative  terras  are  needed,  in  general,  for  a  significant  five- 
digit  accuracy.  The  Importance  of  Equation  (2. 1.3-3)  Is  that  It  makes  the 
sampled-data  frequency  responses  of  the  continuous  portion  of  the  system 
obtainable  without  performing  the  z-transform  operations  (a  significant 
savings  In  computational  effort). 

The  other  frequency  responses  needed  to  complete  the  evaluation  of 
(2. 1.3-2)  at  specified  frequencies  are  of  the  digital  controller  transfer  func¬ 
tion  matrices,  Di(Z)  and  02(Z).  Since  the  elements  of  these  matrices  are 


assumed  to  be  z-transfer  functions,  their  frequency  responses  can  be  evalu¬ 
ated  directly  by  letting  z  =  ej“T,  i.e., 

D{eJ“T)  =  D(z)| 


2  = 


(2. 1.3-4) 


where 


(I)  is  frequency  in  rad/s 
T  is  sample  period  in  seconds 
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e  is  the  natural  log  base. 

Fig.  4  implies  that  both  forward  and  feedback  path  digital  controllers 
are  to  be  used;  however,  in  general  only  one  or  the  other  is  needed.  In  the 
ICAT  design  process  the  specification  of  [DiCz)],  [D2(z)]  and  R(z)  provide 
efficient  means  of  opening  and  closing  loops  and  determining  specific  open 
loop  frequency  responses. 

For  example,  if  [Oi(z)]  =  [0]  then  all  loops  are  open.  Specification  of 
a  controller  for  the  row  i,  column  j  element  of  [DjCz)]  with  the  remainder  of 
the  elements  of  [Di(z)]  set  to  zero  results  in  the  closure  of  a  loop  from 
the  jth  output  to  the  ith  input  through  the  digital  controller  represented  by 
the  specified  element.  Simultaneous  or  independent  closure  of  as  many  loops 
as  desired  may  be  accomplished  in  this  manner. 

In  the  computation  of  open  loop  frequency  responses,  the  input  from 
which  responses  are  being  computed  is  acknowledged  by  specification  of  a 
value  of  I+jO  over  all  frequency  for  the  element  of  R(z)  corresponding  to 
that  particular  input,  with  the  remaining  elements  of  R(z)  set  equal  to  zero. 
The  elements  of  the  sampled-data  output  vector,  C*(s),  are  then  the  fre¬ 
quency  responses  from  the  specified  input  to  each  of  the  outputs  represented 
by  the  elements  of  C*(s).  This  is  to  say  that  each  time  C*{s)  is  computed. 
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a  set  of  tn  frequency  responses  is  generated,  any  or  all  of  which  may  be  of 
interest. 


2.1.4  Application  of  Candidate  Controller  to  ASCOT-I  Structural  Model. 

Introduction 

This  section  describes  the  application  of  the  S2D2  candidate,  ICAT,  to 
the  ASCOT-I  structural  model  as  a  test  case.  A  full  description  of  the 
ASCOT-I  structure  is  given  in  Section  2.2.2  along  with  development  of  the 
dynamic  model.  The  following  paragraphs  provide  a  description  of  the  design 
model  and  the  chosen  control  system  architecture.  A  step  by  step  examination 
of  the  loop  closures  and  resulting  stability  margins  is  then  given  along  with 
a  summary  of  the  compensation  used.  Finally,  conclusions  are  presented. 

ASCOT-I  Design  Model.  The  ASCOT-I  design  model  is  derived  from  the  full 
order  structural  model  of  the  LSS  Ground  Test  Verification  (LSS/GTV)  facility 
structure  housed  at  Marshall  Space  Flight  Center  (MSFC),  Huntsville,  Alabama. 
Use  of  this  model  was  motivated  by  the  desire  to  apply  S2D2  to  a  real ,  exist¬ 
ing  LSS  type  structure.  A  complete  description  of  the  structure  and  accom¬ 
panying  facility  appears  in  Section  2.2.2. 

Modal  truncation  on  the  full  order  dynamic  model  is  accomplished  by 
simply  including  several  modes  with  frequency  less  than  the  half  sample  fre¬ 
quency  of  the  digital  controller  (25  Hz.).  The  simplicity  of  this  structure 
and  the  resulting  model  alleviates  the  need  for  a  more  sophisticated  trun¬ 
cation  technique.  The  resulting  design  model  includes  three  zero  frequency 
rigid  body  modes,  two  "pendulum"  modes,  and  seven  flexible  body  modes.  The 
zero  frequency  modes  represent  translation  of  the  system  at  the  Base 
Excitation  System  (BET)  (two  degrees  of  freedom),  and  rotation  about  the 
z-axis  gimbal .  The  pendulum  modes  are  actually  rigid  body  in  shape  (rotation 


about  X  and  y)  but  hai»«  non-zero  vibrational  frequencies  because  of  the 
effects  of  gravity  upon  the  system.  The  remaining  seven  modes  involve  beam 
bending.  Because  of  the  symmetry  of  the  structure,  this  model  has  no  cross 
axis  coupling.  The  ICAT  technique  can  hand! e  systems  which  are  coupled 
between  axes.  An  example  of  application  of  ICAT  to  such  a  system  is  given  in 
Section  2.2.5. 

The  effector  inputs  for  this  design  are  torques  applied  to  the  structure 
by  the  three  axis  gimbal  system.  The  measurements  are  the  angular  rates  of 
the  Advanced  Gimbal  System  (AGS)  faceplate  about  the  x,  y,  and  z  axes. 

Fig.  5  presents  a  block  diagram  of  the  ASCOT-I  control  system  in  the 
x-axis.  The  y-axis  diagram  is  identical  to  Fig.  5  except  that  each  "x"  is 
replaced  by  "y"  and  vice  versa.  The  z-axis  model  is  similar  also,  except 
that  the  disturbance  force  input  and  line  of  sight  (LOS)  output  are  not 
defined  for  the  z-axis.  The  LOS  measurement  is  computed  as  described  in 
Section  2.2.3.  In  each  axis,  feedback  of  the  angular  rate  of  the  faceplate 
(©*)  is  used  to  add  damping  to  the  structural  vibrations.  The  "outer"  loop 
in  each  axis  uses  derived  faceplate  position  (e*)  to  effect  a  position  command 
loop.  This  type  loop  is  considered  necessary  for  pointing  applications; 
however,  note  that  faceplate  angle  is  the  controlled  quantity  rather  than  LOS 
directly.  It  is  assumed  in  this  analysis  that  LOS  cannot  be  measured 
directly.  The  numerical  integration  used  to  derive  faceplate  position  is 
accomplished  in  the  strapdown  algorithm  (See  Section  2.2.2)  using  quaternions 
such  that  all  measurements  used  for  control  are  with  respect  to  the  labora¬ 
tory  reference  frame.  For  the  purposes  of  control  system  design,  the  numeri¬ 
cal  integration  can  be  represented  by  its  z-transform  and  the  coordinate 
transformations  can  be  ignored. 


STmiCTlIRE 


Design  of  the  ASCOT-I  Controller.  The  design  of  the  compensation  desig¬ 


nated  by  D0  (z)  and  ZaAz)  is  begun  by  examination  of  the  transfer  character- 
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istic  of  ©x/eix  which  is  shown  in  the  Bode  plot  of  Fig.  6.  In  exercising  the 
ICAT  technique,  loops  are  typically  closed  in  order  of  decreasing  control 
bandwidth.  Examination  of  Fig.  6  indicates  a  phase  response  which  lies 
between  ±90“  over  much  of  the  range  of  interest,  which  is  expected  for  a  loop 
using  a  rate  measurement  and  a  collocated  actuator.  This  gives  the  rate  loop 
promise  of  being  the  widest  bandwidth  of  the  two,  which  it  is  indeed.  Also, 
this  response  lacks  classical  "rigid  body"  behavior  even  though  the  sensor 
and  actuator  can  sense  and  effect  this  type  of  behavior.  This  is  because  the 
"rigid  body"  mode  acts  with  a  "pendulum"  behavior,  as  described  earlier,  and 
reveals  itself  as  what  appears  to  be  a  flexible  body  mode  at  about  one 
rad/sec.  The  pendulum  behavior  caused  by  gravity  is  the  thing  which  must 
actually  be  controlled  in  the  position  loop,  and  this  plays  heavily  in  the 
design  and  resulting  performance.  It  is  considered  by  the  author  to  be  a 
worse  case  than  that  in  zero  g. 

In  choosing  compensation  for  the  rate  loop,  it  is  desired  that  as  many 
magnitude  peaks  of  the  compensated  loop  frequency  response  as  possible  peak 
above  OdB.  This  will  introduce  Modal  Suppression  by  Phase  Stabilization 
(MSPS)  [1].  The  higher  a  modal  peak  is  above  OdB,  the  more  will  be  the 
damping  for  that  mode  in  the  transfer  characteristic  of  the  position  loop 
with  the  rate  loop  closed.  Obtaining  MSPS  for  the  flexible  body  modes  is 
easily  accomplished  with  a  simple  feedback  gain  in  D0j^{z);  however,  MSPS  for 
the  pendulum  mode  requires  a  loop  gain  of  sufficient  magnitude  that  dynamic 
compensation  is  required  in  D©  (z).  Damping  of  the  pendulum  mode  is  important 


because,  as  mentioned  before,  this  is  the  motion  to  be  controlled  by  the 
position  loop. 

Using  techniques  described  in  [2],  the  rate  loop  compensation  (Dqj^(z)) 
is  designed  and  is  the  cascade  combination  of  a  first  order  dominant  pole,  a 
lead  compensator,  a  lag  compensator,  and  a  55dB  gain  factor.  The  frequency 
response  of  is  shown  in  Fig.  7,  and  its  mathematical  form  is  summa¬ 

rized  along  with  that  of  the  remainder  of  the  compensation  in  TABLE  1. 

The  compensated  rate  loop  frequency  response  (6x/eix*Dej^(z) )  is  shown  in 
the  Bode  plot  of  Fig.  8  and  in  the  polar  plot  of  Fig.  9.  Application  of  the 
Nyquist  stability  criterion  to  Fig.  9  reveals  a  stable  closed  loop  system 
with  about  15dB  of  gain  margin  and  24°  of  phase  margin. 

The  position  loop  design  begins  with  examination  of  @x/®2x  shown  in 
Fig  10.  Note  the  reduced  height  of  the  modal  peaks  in  this  plot.  An  inte¬ 
gration  is  immediately  specified  as  part  of  the  position  loop  compensation 
{E0jj(z))  so  that  the  system  will  be  Type  I  and  therefore  have  no  steady  state 
error  due  to  a  constant  command  input.  Inclusion  of  the  integrator  charac¬ 
teristic  in  the  position  loop  frequency  response  results  in  the  Bode  plot  of 
Fig.  11.  The  loop  could  be  stabilized  at  this  point  using  a  gain  factor  in 
addition  to  the  integrator  already  chosen;  however,  the  bandwidth  would  be 
very  limited  and/or  the  stability  margins  small.  For  this  reason,  a  second 
order  lead  compensator  having  an  underdamped  zero  pair  and  an  underdamped 
pole  pair  is  designed  to  phase  stabilize  the  pendulum  mode  and,  as  a  result, 
enhance  the  bandwidth  of  the  position  loop  without  sacrificing  stability 
margins.  The  frequency  response  of  the  resulting  compensation  is  shown  in 
Fig.  12.  The  frequency  range  over  which  the  compensation  exhibits  the  large 
phase  lead  required  to  phase  stabilize  the  pendulum  mode  is  narrow,  and  use 
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Figure  10.  Frequency  Response  of  Open  X-Axis 
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Figure  11.  Frequency  Response  of  Open  X-Axis  Position  Looi 
Cascaded  with  an  Integrator. 


of  such  compensation  vrould  typically  be  avoided  because  of  robustness  problems 
resulting  from  uncertainty  of  the  modal  frequency.  However,  the  pendulum 
frequency  is  quite  predictable  and  very  consistent;  therefore,  compensation 
such  as  this  is  acceptable  in  this  situation. 

The  position  loop  frequency  response  with  the  rate  loop  closed  and  all 
compensation  included  (both  rate  loop  and  position  loop)  is  shown  in  the  Bode 
plot  of  Fig.  13  and  the  polar  plot  of  Fig.  14.  Examination  of  Fig.  14  indi¬ 
cates  a  gain  margin  of  lOdB  and  a  phase  margin  of  64‘. 

Design  of  compensation  for  the  y-axis  control  system  is  similar  to  that 
of  the  x-axis,  and,  in  fact,  the  compensation  differs  only  by  the  gain  fac¬ 
tors  used.  The  y-axis  compensation  along  with  the  resulting  stability 
margins  is  summarized  in  TABLE  1. 

The  z-axis  control  system  has  an  architecture  similar  to  that  of  the  x 
and  y-axes  and  is  designed  in  much  the  same  manner.  However,  it  is  of  much 
less  interest  here  for  three  reasons:  (1)  it  is  not  a  pointing  axis,  (2)  it 
in  no  way  affects  the  LOS  pointing  because  of  the  decoupled  nature  of  the 
dynamic  model,  and  (3)  only  one  flexible  body  mode  appears  in  the  frequency 
range  of  interest,  and  therefore  the  design  is  relatively  simple.  For  these 
reasons,  the  z-axis  design  is  not  discussed  in  detail;  however,  the  compensa¬ 
tion  and  resulting  stability  margins  are  included  in  TABLE  1. 

Conclusions  regarding  the  ASCOT-I  design.  Design  of  a  three  axis  con¬ 
trol  system  for  the  ASCOT-I  model  using  the  candidate  S2D2,  ICAT,  has  been 
demonstrated.  The  resulting  design  produces  a  stable  closed  loop  system  with 
a  reasonable  amount  of  robustness  as  indicated  by  the  stability  margins  of 
TABLE  1.  The  ASCOT-I  model  does  not  provide  a  very  challenging  problem  for 
ICAT  because  of  the  decoupled  nature  of  the  model.  It  does,  however,  provide 


the  reader  with  an  introduction  to  the  step  by  step  nature  of  ICAT  and  the 
accompanying  stability  checks.  A  more  thorough  application  of  ICAT  is  pro¬ 
vided  in  Section  2.2.5.  Section  2.1.6  discusses  the  performance  and  time 
domain  stability  of  the  ASCOT-I  design. 

2.1.5  Design  of  Prefilters. 

Introduction.  The  effect  of  unmodeled  high  frequency  system  components 
is  often  neglected  in  the  design  of  sarapled-data  control  systems.  The 
rationale  for  leaving  these  components  unmodeled  is  that  they  will  either 
not  be  excited  or  will  contribute  very  little  to  the  low  frequency  charac¬ 
teristics. 

The  advent  of  Large  Space  Structures  (LSS)  with  significant  high  fre¬ 
quency  flexible  modes  introduces  almost  a  worst  case  of  unmodeled  high  fre¬ 
quency  components.  A  legitimate  concern  for  control  systems  designers  for 
such  structures  is  whether  neglecting  the  high  frequency  modes  contributes  to 
system  performance  degradation. 

Since  control  systems  for  LSS  are  envisioned  to  be  sampled-data  in 
nature,  the  main  problem  expected  to  occur  is  that  of  aliasing.  The  phenome¬ 
non  of  aliasing  is  most  easily  understood  by  considering  the  Laplace 
Transform  of  the  output  of  an  ideal  sampler  with  an  input  signal  whose 
Laplace  Transform  is  X{s)  (See  Fig.  15.). 

It  can  be  shown  that 
1  - 

X*(s)  =  -  I  X(s  +  jna,s) 

T  n=—  (2. 1.5-1) 


where  uj  is  the  sampling  frequency  in  rad/sec.  Figure  16  is  a  representation 
of  the  s-plane  with  the  poles  of  X{s)  shown  schematically. 


lii 


Examination  of  Equation  (2. 1.5-1)  results  in  the  conclusion  that,  if  X(s) 
has  poles  at  oj  ^  j“s/4  ®2  i  J“s3/4  then  X*(s)  has  poles  at  ai  jusM  + 

jnuis  and  02  ^  j(i)s3/4  ^  jna)$.  jp  particular  X*{s)  has  poles  at  02  J  “s/4 
and  j  ajs/4,  the  significance  being  that  two  poles  are  now  within  the 

primary  strip  (-015/2  <  oi  <  015/2)  that  do  not  appear  there  for  the  unsampled 
X(s).  These  two  poles  are  said  to  have  aliased  back  into  the  primary  strip 
and  can  cause  serious  stability  problems  since  a  control  system  containing  an 
ideal  sampler  may  introduce  low-frequency  poles  that  the  designer  has  not 
considered. 

The  problem  of  aliasing  is  usually  solved  in  two  ways.  First  the 
designer  may  know  that  his  system  is  effectively  bandlimited  and  choose  the 
sampling  period  so  that  015/2  is  greater  than  the  highest  significant  fre¬ 
quency,  Alternatively  the  designer  may  decide  that  the  additional  hardware 
required  to  implement  a  very  fast  sampled-data  control  system  is  not  justi¬ 
fied.  In  this  case  analog  prefiltering  is  used  to  bandlimit  the  system  to  a 
reasonable  frequency  range  by  inserting  a  low  pass  filter  with  a  cutoff  fre¬ 
quency  chosen  so  that  the  undesirable  high  frequency  components  are  signi¬ 
ficantly  attenuated  at  frequencies  higher  than  015/2. 

Two  alternatives  to  the  approach  of  analog  prefiltering  are  presented 
here.  The  first  method  involves  using  a  digital  prefilter  with  a  faster  uni¬ 
form  sampling  rate.  The  second  method  involves  a  non-uniform  sampling  rate 
with  an  average  rate  slower  than  that  used  by  a  corresponding  fast  uniform 
rate  prefilter. 

Different  justifications  exist  for  each  of  the  two  approaches.  In  the 
case  of  the  fast  uniform  rate  prefilter,  the  additional  hardware  cost  can  be 
kept  comparitively  low  by  the  use  of  dedicated  hardware. 


The  non-uni  form  rate  prefilter,  if  proven  practicable,  has  the  advantage 
that  extremely  fast  sampling  is  not  required. 

Before  examining  the  two  prefilter  schemes  in  detail,  it  is  worthwhile 
to  consider  the  circumstances  in  which  aliasing  may  lead  to  stability  and/or 
performance  degradation.  The  most  efficient  method  of  examining  this  problem 
is  by  considering  the  s-plane  representation  of  pole  locations. 

Consider  the  basic  sampled-data  feedback  system  of  Figure  17.  Suppose 
G(s)  has  2  poles  inside  the  primary  strip  and  2  poles  outside  the  primary 
strip.  The  poles  of  G{s)  are  shown  on  the  s-plane  in  Figure  18.  Compensator 
design  is  generally  effected  in  a  sampled  data  system  through  use  of  the 
sampled  data  loop  frequency  response.  In  this  case  the  loop  response  is 
G*(s),  whose  poles  are  represented  in  Figure  19.  Two  poles  have  now 
“aliased  back"  into  the  primary  strip  where  none  had  previously  appeared. 
Now  the  system  stability  is  determined  by  the  sampled-data  characteristic 
equation: 

A(s)  =  1  +  G*(s)  =  0  (2. 1.5-2) 

Clearly,  G*{s)  is  periodic  with  period  juj.  Hence  G*(ju))  is  completely 
described  by  the  portion  of  the  ju  axis  between  -0)5/2  and  +0)5/2.  Normally, 
the  periodicity  of  the  poles  of  G*{jo))  contribute  very  little  to  its  value  in 
this  range,  as  G(jo))  is  assumed  to  have  no  significant  values  for  frequencies 
above  0)5.  However,  if  this  is  not  true,  which  is  the  case  when  aliasing 
occurs,  G*(jo))  can  be  affected  so  as  to  cause  a  degradation  in  stability.  It 
should  be  pointed  out  that  an  actual  aliasing  of  poles  is  not  required  for  a 
stability  degradation.  In  fact,  whenever  G(jo))  has  significant  components  at 
frequencies  higher  than  0)5/2  stability  can  be  affected.  It  is  for  this 


reason  that  the  effect  of  decreasing  stability  margins  in  a  sampled-data 
system  is  termed  "frequency  folding." 

A  topic  of  particular  interest  to  designers  of  Large  Space  Structures  is 
that  of  lightly  damped  modes  which  may  unavoidably  occur  at  frequencies 
higher  than  uj/Z.  While  the  existence  of  such  modes  has  been  the  subject  of 
great  concern,  it  should  be  noted  that  when  these  modes  are  aliased  below 
(i>s/2,  there  exists  a  bound  on  the  damping  ratio,  based  on  the  fact  that  the 
real  part  of  the  resulting  aliased  poles  cannot  change  whereas  the  imaginary 
part  of  the  poles  may  be  much  reduced. 

Fast  Uniform  Rate  Prefilters.  The  effect  of  prefilters  on  aliasing  phe¬ 
nomenon  can  be  considered  from  either  a  signal  or  loop  gain  point  of  view. 
In  the  signal  point  of  view,  the  designer  inserts  a  prefilter  in  order  to 
eliminate  the  excitation  of  undesirable  high  frequency  system  characteristics. 
For  the  loop  gain,  or  feedback,  point  of  view  the  designer  considers  the 
function  of  the  prefilter  as  reducing  the  loop  bandwidth  of  the  system,  hence 
eliminating  the  high  frequency  system  components  of  concern.  Both  of  the 
above  mentioned  ways  of  visualizing  the  function  of  an  anti-aliasing  pre¬ 
filter  are  equally  valid  and  both  concepts  can  be  used  interchangeably. 

The  concept  of  digital  prefiltering  arises  naturally  in  the  scenario  of 
on-orbit  control  system  redesign.  For  example,  earth  bound  designers  may 
desire,  after  in-flight  testing,  to  incorporate  in  the  control  system  design 
previously  unmodeled  modes.  These  modes,  if  they  occur  at  frequencies  above 
the  cut-off  frequency  of  the  prefilters,  can  be  uncontrollable.  An  analog 
prefilter  with  its  fixed  cutoff  frequency  would  not  allow  the  designers  to 
take  advantage  of  the  in-flight  data.  On  the  other  hand,  designers  can 
easily  change  the  cutoff  frequency  of  a  digital  prefilter  from  the  ground. 
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However,  from  the  previous  discussion  of  the  effect  of  sampling  it  should  be 
clear  that  a  digital  prefilter  cannot  filter  at  frequencies  above  015/2,  a 
direct  consequence  of  the  periodicity  of  the  sampled  data  frequency  response. 

The  unfortunate  fact  that  the  cutoff  frequency  of  uniform  rate  digital 
prefilters  is  limited  by  015/2  leads  swiftly  to  the  conclusion  that  the 
sampling  rate  associated  with  the  prefilter  roust  be  very  high  to  ensure  that 
all  significant  high  frequency  components  are  attenuated.  The  practicality 
of  such  a  filter  is  determined  by  whether  it  is  feasible  to  implement  such  a 
high  rate  filter  with  available  hardware  and  whether  it  can  be  assumed  that 
the  system  can  be  considered  effectively  bandlimited  to  a  reasonable  fre¬ 
quency  range. 

If  the  sampling  rate  is  limited  only  by  the  speed  at  which  digital 
filter  calculations  can  be  made  then  the  possibility  exists  that  a  dedicated 
hardware  digital  prefilter  is  practical.  On  the  other  hand,  if  the  sampling 
rate  is  limited  by  the  availability  of  data,  the  probability  is  that  a  fast 
uniform  rate  digital  prefilter  will  be  impractical. 

The  design  of  a  fast  uniform  rate  digital  prefilter  is  straightforward. 
First,  the  fast  sampling  rate  associated  with  the  digital  prefilter  should, 
if  possible,  be  at  least  ten  times  the  sampling  rate  associated  with  the 
slower  digital  controller  in  order  to  minimize  the  fast  rate  sampler's  effect 
on  the  overall  system  design.*  Second,  the  cut-off  frequency  of  the  prefilter 
should  be  selected  to  be  at  or  below  the  half  sample  rate  of  the  controller. 
Commonly  a  simple  second  order  low  pass  characteristic  is  chosen,  with  its 


*lf  the  fast  rate  sampler  is  too  slow,  then  multi  rate  design  techniques 
may  be  required  [3]. 
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pair  of  complex  conjugate  poles  critically  damped  at  a  frequency  of  (us/2 
(slow  rate).  Such  a  filter  offers  an  acceptable  tradeoff  between  gain  and 
phase  contribution,  translating  to  a  tradeoff  between  system  stability  and 
system  response  time. 

Non-Uniform  Rate  Sampling.  The  ideal  marriage  of  the  fast  rate  digital 
prefilter  and  the  analog  prefilter  is  a  prefilter  which  is  digital  in  nature 
while  at  the  same  time  not  requiring  a  very  fast  sampling  scheme.  Such  a 
filter,  while  not  possible  with  standard  uniform  rate  sampling  schemes,  may 
be  possible  if  unconventional  sampling  schemes  are  used.  In  order  to  see 
that  aliasing  is  closely  related  to  the  sampling  scheme,  the  sampling  process 
is  modeled  as  an  impulse  modulator,  as  in  Figure  20,  with  p(t)  a  periodic 


train  of  impulses,  i.e.. 


p(t)  =  I  6(t  -  nT). 
n=0 


It  is  well  known  that 


(2. 1.5-3) 


i{x(t)p(t)}  *  X(s)*P(s)  (2. 1.5-4) 

where  J[.{  }  denotes  the  Laplace  Transform  operator,  *  denotes  the  operation 
of  complex  convolution,  and  X(s)  and  P(s)  are  the  Laplace  transforms  of  x(t) 
and  p(t),  respectively.  The  Laplace  transform  of  p(t)  can  easily  be  evalu¬ 


ated  as 


P(s)  =  I  e-nsT  =  - 

n=0  1  -  e-sT  . 


(2. 1.5-5) 


From  Equation  4 


go 

1  f  1 

X*(s) - /  X{u) - - - - 

-'c-j-  1  -  e-T(s-u) 


I  residues  of 


1  _  e-T(s-u) 


(2. 1.5-6) 
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Now  there  are  poles  of  - 

1  .  e'T(u-s) 


at  u  =  s  ±  jnus,  n  =  0,  1, 


•  •  • 


and  the 


corresponding  residues  are  X(s  ±  jnuisl/T. 
1  • 

X*(s)  =  -  I  X(s  +  jnus). 

T  n=-- 


Hence  X^(s)  can  be  written 


(2. 1.5-7) 


and,  as  discussed  previously,  the  aliased  poles  and  frequency  folding  occur 
because  of  the  sidebands  generated  by  Equation  (2. 1.5-7).  It  should  be  clear 
from  the  development  of  Equation  (2. 1.5-7)  that  these  effects  are  immediate 
consequences  of  the  periodicity  of  the  poles  of  P(s),  which  in  turn  results 
from  the  periodicity  of  the  sampling  impulse  train  p(t). 

Because  the  sidebands  which  cause  aliasing  and  frequency  folding  are 
results  of  the  periodicity  of  the  sampling  process,  an  alternate  non-periodic 
sampling  scheme  is  proposed  which,  at  least  theoretically,  does  not  cause 
aliasing.  This  alternate  sampling  scheme  has  the  desirable  feature  of  having 
a  minimum  time  between  samples,  thus  allowing  implementation  without  the 
unrealistic  requirement  of  obtaining  samples  arbitrarily  close  together  in 
time.  A  digital  prefilter  is  included  to  investigate  the  possibility  of 
filtering  above  the  "average"  half  sample  rate  associated  with  the  sampling 
process. 

The  non-uniform  rate  sampling  process  is  defined  by  the  following  steps. 

1.  A  sample  is  assumed  to  be  obtained  at  t  -  0. 

2.  The  next  sample  time  is  determined  from  a  random  number  generator 


uniformly  distributed  on  the  interval  [Tnjin.  Tnjax3* 
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The  resulting  sample  of  the  random  number  generator  is  added  to  the  pre¬ 
vious  sample  time  to  determine  the  time  of  the  next  sample.  Mathematically 
the  sampling  process  is  again  modeled  as  an  impulse  modulator  with  impulse 
train  p{t)  defined  by 


p(t)  =  6[t  -  (T  +  ...  +!„)]  (2. 1.5. 8) 

where  the  T-j's  are  independent  random  variables  uniformly  distributed  on  the 
interval  [Tmin,  Tmax]- 

In  order  to  prove  that  the  sampling  process  does  not  introduce  side¬ 
bands,  an  expression  for  the  output  power  spectral  density  (PSD)  in  terms  of 
the  input  PSD  must  be  derived.  Since  the  output  of  the  sampler  is  the  pro¬ 
duct  of  the  input  x{t)  and  the  modulating  impulse  train  p(t)  the  autocorrela¬ 
tion  function  of  the  output  y(t)  can  be  written  as  follows: 
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Ryy(T)  =  E[y(t  +  T)y(t)]  »  E[x(t  +  T)p(t  +  T)x(t)p(t)]  (2. 1.5-9) 


Since  x(t)  and  p(t)  are  statistically  independent, 

Ryy{T)  =  ECx(t  +  T)x{t)]  E  [p{t  +  T)p(t)]  =  Rxx(t)  Rpp(y)  (2.1.5-10) 
Then  the  PSD  of  the  output  is 


Sy(a))  =  Sx((u)*Sp(a»)  (2.1.5-11) 

since  the  PSD  is  simply  the  Fourier  transform  of  the  autocorrelation  func¬ 
tion.  An  important  fact  to  note  is  that  if  Sp(a))  =  6((d),  then  Sy(u)  =  Sx(u)); 
a  very  desirable  result  since  there  are  no  sidebands  in  this  case. 

Equation  11  implies  that  if  Sp((jo)  is  a  good  approximation  of  6((d),  then 
aliasing  should  not  be  a  significant  probl«n.  The  PSD  of  p(t)  can  be  found 
by  observing  that  [4] 
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where  T  Is  the  time  interval  of  p(t)  used  to  obtain  the  Fourier  transform 
a  given  T,  there  is  always  a  k  such  that 
Px(t)  =  6(t)  +  6(t  -  Ti)  +  ...  +  6(t  -  Ti  -  T2  ...  -  T^)  (2.1.5.13) 

so  that  the  Fourier  transform  of  Pj(t)  is 

PjCju)  =  1  +  e-j“Ti  +  ...  +  +  ...  T^.i).  (2.1.5.14) 


Hence, 
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2.1.5-15 


Equation  (2.1.5-15)  can  be  written  as 


-  PTCj*-) 
T 


Tj  +  ...  +  T|j 


Tj  +  ...  +  T|j 


+  ...  +  T|( 


[eJoiTj  +  +  eJudi  +  ...  +  Tk.i)] 


[e4<»T2  +  +  ejw(T2  +  ...  T|c-i)] 


+  ...  +T|j 


+  ...  +  T|j 


ej“T|(_i 


[e-jwTi  +  ...  +  +  ...  +  T^.i)] 


e-juTij.i  . 


Ti  +  ...  +  Tv 
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Now,  as  T  gets  large - approaches -  ,  where  is  the  mean 

Ti  +  ...  +  T|^  kT^ 

value  of  each  of  the  T^-'s.  The  required  expectation  can  be  found  terra  by 
term  from  Equation  (2.1.5-16).  Since  the  T^'s  are  identically  distributed, 

ECejuTf]  =  E[ej“Tj]  f^ir  any  i  and  j.  (2.1.5-17) 


Also,  since  the  T^'s  are  independent, 

£[ej<»»(Ti  +  ...  +  T|(_i)]  =  ECej“*^l^  •••  E[ej<**^k-1^ 

=  E[ej“Ti]k-l. 

The  probability  density  function  of  each  T^  is 

,  ,  .  fl/f^max  ~  ^min^  J  ^min  ^i  i.  ^max 
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;  otherwise. 
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(2.1.5-18) 
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(2.1.5-20) 


Xn  =  E[e-J'*'Ti] 
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Substituting  into  Equation  (2.1.5-16)  and  gathering  terms. 
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(2.1.5.21) 


(2.1.5-22) 


Taking  the  limit  as  k  »  gives 
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(2.1.5-23) 


Investigation  of  Equation  (2.1.5-23)  yields  the  desirable  result  that  Sp((u) 
is  indeed  an  approximation  of  an  impulse  centered  at  zero  with  a  constant 


is  expected  that  the  non-uniform  rate  sampler  greatly  reduces  aliasing  phe¬ 
nomenon.  However,  the  question  of  whether  a  suitable  digital  prefilter  is 
possible  using  the  sampling  scheme  is  unresolved.  A  theoretical  discussion 
of  the  digital  prefilter  using  a  non-uniform  rate  sampling  scheme  is  very 
difficult,  if  not  impossible,  and  is  not  Included  here  except  to  suggest 
future  study  on  the  matter.  However,  computer  simulations  with  such  a  pre¬ 
filter  were  carried  out  in  order  to  assess  its  practicality.  The  system 
simulated  along  with  the  simulation  results  are  presented  in  the  following 
sections. 

Non-Uniform  Rate  Digital  Prefilter  Simulation  Results.  Computer  simul a- 
tion  results  for  a  typical  system  subject  to  aliasing  problems  are  included. 
The  non-uniform  rate  digital  prefilter  is  observed  to  improve  the  system  per¬ 
formance  only  as  concerns  the  effects  of  aliasing.  The  overall  system 
response  exhibits  a  degradation  in  stability  margins. 

Assume  that  one  is  required  to  control  a  double  integrator  system  with  a 
transfer  function  given  by 
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term  l/T^.  From  the  previous  discussion  of  the  desirable  form  of  Sp(u),  it 
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(2.1.5-24) 


Also,  assume  that  this  particular  system  exhibits  an  unknown  vibrational  mode 
so  that  the  actual  plant  transfer  function  is  given  by 


i  w  i 


. -I 

••  -■  ■.< 


'.^1 


• '  •  ■  •  •  .  I 


s  s  +  2i:ain5  +  (2.1.5-25) 

A  block  diagram  for  this  system  is  shown  in  Fig.  21.  In  this  block  diagram 

Di(z)  is  a  digital  compensator  with  a  uniform  rate  sampler  and  D2{z)  is  a 
digital  prefilter  with  a  non-uniform  rate  sampler.  Note  that  the  compensator 
precedes  the  prefilter  in  the  forward  path  of  this  system;  however,  this  need 
not  be  the  case. 

Since  the  vibrational  characteristics  of  this  plant  are  unknown  to  the 
system  designer.  Equation  (2.1.5-24)  will  be  used  to  choose  a  suitable  uniform 

sample  rate  and  design  the  compensator  Di(z).  In  this  example  the  uniform 

sample  rate  is  chosen  to  be  200  rad/sec.  The  gain  kj  in  Equation  24  is 
chosen  to  be  177.8  so  that  the  sampled-data  frequency  response  of  the  plant 
crosses  zero  dB  at  approximately  13  rad/sec.  A  lag-lead  compensator  is  now 
designed  so  that  the  loop  gain  crossover  frequency  remains  13  rad/sec  with  a 
phase  margin  of  40  degrees  and  a  gain  margin  of  about  12  dB.  The  transfer 
function  for  this  compensator  is  given  by 

.3019  z  -  .2864  7.039  z  -  6.21 

Oi(z)  =  -  -  . 

z  -  .9845  z  -  .1707  (2.1.5-26) 

Fig.  22  shows  a  Bode  plot  of  the  compensated  loop  frequency  response  and  Fig. 

23  shows  a  unit  step  response  for  the  compensated  system. 

For  the  purposes  of  the  example  it  is  now  necessary  to  choose  the  dc 
gain,  damping  ratio,  and  natural  frequency  for  the  unmodeled  vibrational  mode 
so  that  an  alias  occurs  in  such  a  way  as  to  adversely  effect  the  stability  of 
the  system.  Two  such  cases  are  chosen  for  this  example.  In  both  cases  the 
alias  frequency  is  50  rad/sec.  For  Case  1  the  natural  frequency  is  chosen  to 
be  150  rad/ sec,  the  damping  ratio  is  .01,  and  the  gain  k?  is  .02.  A  Bode 
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plot  of  the  sampled-data  loop  frequency  response  of  the  compensated  system 
with  this  vibrational  mode  is  shown  in  Fig.  24.  The  resonant  peak  due  to  the 
alias  results  in  a  very  low  (approximately  3  dB)  gain  margin  so  that  the 
closed  loop  system  should  exhibit  an  oscillatory  time  response. 

For  Case  2  the  vibrational  mode  natural  frequency  is  chosen  to  be  350 
rad/ sec,  the  damping  ratio  is  .005,  and  k2  is  .02.  A  Bode  plot  of  the  loop 
frequency  response  of  the  compensated  system  with  this  vibrational  mode  is 
shown  in  Fig.  25. 

Fig.  26  shows  the  unit  step  response  of  the  system  with  the  150  rad/sec 
vibrational  mode,  and  Figure  27  shows  the  response  of  the  system  with  the  350 
rad/ sec  vibrational  mode.  The  oscillatory  behavior  of  both  of  these  respon¬ 
ses  illustrates  the  need  for  a  prefilter  to  reduce  the  effects  of  the  alias. 

The  prefilter  should  be  designed  in  such  a  way  as  to  preserve  the  sta¬ 
bility  margins  obtained  with  the  compensator  designed  above.  For  this 
example  a  second  order  low  pass  filter  with  a  natural  frequency  of  88  rad/sec 
and  a  damping  ratio  of  .3  is  chosen.  This  filter  should  contribute  less  than 
5  degrees  of  phase  lag  at  the  loop  gain  crossover  frequency  (13  rad/sec)  and 
only  about  3  dB  of  gain  at  the  loop  phase  crossover  frequency  (40  rad/sec). 
Therefore,  the  loop  stability  margins  are  not  greatly  effected  by  this 
filter.  Also,  the  filter  should  provide  about  8  dB  of  attenuation  at  150 
rad/sec  and  about  25  dB  at  350  rad/sec.  This  attenuation  should  greatly 
reduce  the  oscillatory  effects  of  the  alias.  The  transfer  function  for  the 
analog  pre-filter  is  given  by 


The  digital  pre-filter  will  be  designed  in  the  w  plane  using  an  average 
sample  period  equal  to  the  uniform  sample  period  associated  with  the  digital 
compensator  (.0314  sec).  In  the  w  plane  the  transfer  function  for  the  pre¬ 


filter  is  given  by 


where 


w  4-  2c  (tin  w  +  un 


(tinz  T 

(tin  =  Tan  ( - ) 


(2.1.5-28) 


(2.1.5-29) 


Using  the  natural  frequency  of  88  rad/sec  and  the  damping  ratio  of  .3  as  in 
the  analog  prefilter  and  using  the  transformation 

z  -  1 
w  = - 

z  +  1  (2.1.5-30) 

yields  the  z  domain  transfer  function  for  the  prefilter  as  shown  below. 

27.885  (z^  +  2z  +  1) 

02(2)  = - 5 - 

32.053  2^  +  53.77  z  +  25.716  (2.1.5-31) 

The  following  section  contains  the  unit  step  responses  of  the  two  sys¬ 
tems  with  both  the  analog  and  the  digital  prefilters  in  the  loop  separately 
The  simulation  was  used  extensively  to  experiment  with  both  the  random  sample 
interval  (T^ax  "  ^min^  placement  of  the  digital  prefilter  in  the 

loop. 

Case  1:  150  rad/ sec  vibrational  mode 

Fig.  28  shows  the  response  of  the  system  with  the  150  rad/sec  vibra¬ 
tional  mode  with  the  analog  pre-filter  in  the  loop.  A  comparison  of  this 
response  to  that  shown  in  Fig.  26  illustrates  the  desirable  effect  of  the 
analog  prefilter  with  regard  to  reducing  the  effect  of  the  alias. 


Figs.  29,  30  and  31  show  the  response  of  the  system  with  the  digital 
prefilter  in  the  loop  as  shown  in  Fig.  21.  The  average  sample  period  for 
the  non-uniform  rate  sampler  is  the  same  for  all  three  responses  and  is  equal 
to  the  uniform  sample  period.  For  the  response  shown  in  Fig.  29  the  random 
sample  rate  was  varied  between  190  rad/sec  and  210  rad/ sec.  The  random 
sample  rate  was  varied  between  150  rad/ sec  and  250  rad/sec  for  the  response 
shown  in  Fig.  30.  For  the  response  shown  in  Fig.  31  the  random  sample  rate 
was  varied  between  100  rad/sec  and  300  rad/sec.  These  responses  show  a 
slight  improvement  in  the  system  response  as  compared  to  the  response  with  no 
prefilter  if  only  the  aliasing  effects  are  considered.  The  150  rad/sec 
oscillations  are  of  smaller  amplitude  and  settle  faster  in  Figs.  29  and  30. 
However,  the  non-uniform  rate  sampling  tends  to  degrade  the  system  stability 
margins.  This  can  be  seen  in  all  three  responses  by  the  lower  frequency 
(approximately  16  rad/sec)  oscillations  which  weren't  present  at  all  with  the 
analog  prefilter.  These  responses  also  show  that  the  system  stability  margins 
get  progressively  worse  as  the  range  of  the  ramdom  sampling  intervals  is 
increased. 

In  the  responses  shown  in  Figs.  32,  33,  and  34,  the  sample  rate  was 
varied  exactly  as  in  Figs.  29,  30,  and  31,  respectively.  The  obvious  differ¬ 
ences  in  these  two  sets  of  responses  are  a  result  of  changing  where  the  non- 
uniform  sampling  rate  prefilter  is  placed  in  the  loop.  The  latter  set  of 
responses  was  generated  with  the  prefilter  placed  before  the  compensator  so 
that  the  random  sampler  is  sampling  the  error  signal.  The  responses  in 
Figs.  32,  33,  and  34  exhibit  basically  the  same  characteristics  as  their 
counterparts  in  Figs.  29,  30,  and  31.  Both  prefilter  configurations  seem 
to  reduce  the  alias  effects  as  well  as  degrade  the  loop  stability  margins. 
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The  latter  configuration  tends  to  produce  step  responses  with  greater  over¬ 
shoot  and  faster  rise  time  than  their  counterparts  in  the  configuration  of 
Fig.  21. 

Case  2:  350  rad/sec  vibrational  mode 

Fig.  35  shows  the  unit  step  response  of  the  system  with  the  350  rad/sec 
vibrational  mode  with  the  analog  prefilter  in  the  loop.  The  analog  prefilter 
is  seen  to  be  more  effective  in  reducing  the  aliasing  effects  for  this  case 
because  the  natural  frequency  is  higher.  The  stability  problems  due  to  the 
alias  practically  disappear  with  the  analog  prefilter  in  the  loop. 

Figs.  36  and  37  show  responses  of  the  system  with  the  non-uniform  samp¬ 
ling  rate  prefilter  in  the  loop  as  shown  in  Fig.  21.  In  the  response  shown 
in  Fig.  36  the  sample  rate  was  varied  between  190  rad/sec  and  210  rad/sec. 
The  sample  rate  was  varied  between  150  rad/ sec  and  250  rad/ sec  in  the  response 
shown  in  Fig.  21.  The  average  sample  rate  for  the  prefilter  is,  therefore, 
the  same  as  the  uniform  sample  rate  for  these  responses.  These  responses 
exhibit  most  of  the  same  characteristics  that  were  present  in  the  responses 
for  Case  1.  Again,  it  is  evident  that  the  prefilter  tends  to  degrade  the 
system  stability  margins,  an  effect  that  worsens  with  increasing  range  of  the 
sampling  intervals.  The  system  responses  with  the  prefilter  moved  ahead  of 
the  compensator  follow  the  same  basic  pattern  as  for  those  shown  for  Case  1 
and  are,  therefore,  omitted. 

Conclusions.  Two  alternatives  to  analog  prefiltering  were  presented  in 
the  preceding  sections.  A  comparison  of  their  respective  advantages  and  dis¬ 
advantages  is  included  in  the  following  discussion.  Future  research  direc¬ 
tions  are  also  recommended. 


with  Bending  Mode  at  350  rad/sec 
log  Prefilter. 
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The  first  of  the  alternatives  to  analog  prefiltering  is  fast  uniform 
rate  digital  prefiltering.  While  no  theoretical  stumbling  blocks  to  the 
implementation  of  such  a  filter  exist,  there  are  two  practical  considerations 
with  which  the  system  designer  must  contend.  The  first  consideration  is  that 
of  whether  appropriate  fast  samples  are  available.  Second,  the  designer 
should  realize  that  the  implementation  of  the  fast  rate  prefilter  will  pro¬ 
bably  entail  the  use  of  dedicated  hardware  and,  hence,  a  higher  cost. 

The  second  alternative  to  analog  prefiltering  is  the  non-uniform  rate 
prefilter.  While  no  exact  theory  exists  to  provide  for  such  a  filter  design, 
the  results  of  Section  III  at  least  show  that  the  non-uniform  rate  sampling 
process  minimizes  the  effects  of  aliasing.  However,  in  Section  lY  it  was 
shown  that  the  use  of  a  corresponding  digital  prefilter  based  on  an  average 
sample  rate  degrades  the  system  stability  margins.  In  spite  of  these  prob¬ 
lems,  future  research  into  the  prefilter  design  is  indicated.  Several  possi¬ 
bilities  for  research  exist.  First,  it  may  be  possible  to  design  compensation 
with  the  non-uniform  rate  sampler,  precluding  the  need  for  a  prefilter  since 
aliases  are  not  introduced.  Second,  a  different  non-uniform  rate  sampling 
scheme  might  yield  better  prefilter  performance.  Indeed,  Leneman's  work  with 
"jitter  rate"  sampling  [5]  appears  especially  promising. 

2.1.6  Investigation  of  System  Stability  and  Performance 


This  section  includes  discussion  of  the  system  performance  of  the 
ASCOT-I  controller  coupled  with  the  ASCOT-I  model  as  described  in  Section 
2.1.4.  The  performance  is  analyzed  with  respect  to  position  loop  bandwidth, 
random  disturbance  rejection,  and  deterministic  disturbance  rejection.  The 
disturbance  rejection  studies  use  the  line  of  sight  (LOS)  criteria  developed 
under  Section  2.2.3.  The  disturbance  inputs  are  described  in  Section  2.1.2. 
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Position  Loop  Bandwidth.  The  position  loop  bandwidth  is  of  interest  as 
an  indication  of  system  performance  in  terms  of  response  to  pointing  command 
inputs  (such  as  a  slew  command).  Figs.  38,  39,  and  40  display  the  dosed 
loop  frequency  responses  from  commanded  position  input  to  measured  faceplate 
position  output  for  the  x,  y,  and  z-axes,  respectively.  The  x-axis  bandwidth 
is  about  0.8  rad/sec  and  the  y-axis  bandwidth  is  about  0.9  rad/ sec.  The 
y-axis  bandwidth  is  wider  because  of  inherent  differences  in  the  gimbal 
positions;  however,  it  exhibits  greater  peaking  near  the  bandwidth.  The 
z-axis  response  of  Fig.  40  is  the  most  well  behaved  of  the  three;  indicative 
of  the  benign  nature  of  the  torsional  behavior.  The  z-axis  bandwidth  is 
roughly  3.0  rad/sec. 

Random  Noise  Disturbance  Suppression.  The  random  noise  suppression 
results  herein  were  generated  with  the  frequency  domain  model  of  the  struc¬ 
ture  which  was  used  in  conducting  the  design.  The  Power  Spectral  Density 
(PSO)  of  the  random  force  disturbance  applied  to  the  structure  is  given  in 
Section  2.1.2.  With  the  coordinate  system  and  input/output  definitions 
chosen  disturbance  force  in  the  y-direction  (Fy)  generates  error  in  x-axis 
LOS  (LOSx)  and  vice  versa.  Cross- axis  disturbance  is  precluded  by  the 
decoupled  nature  of  the  system. 

The  PSD  of  the  LOS  error  in  the  x-axis  with  disturbance  applied  in  the 
y-direction  and  all  control  loops  open  is  given  in  Fig.  41.  Fig.  42  displays 
the  same  quantity  with  all  control  loops  closed.  Comparison  of  these  plots 
shows  a  definite  improvement  in  disturbance  rejection  with  the  control  system 
included,  especially  at  the  modal  peaks.  Quantitatively,  the  RMS  error  in 
the  open  loop  case  is  1.5  X  10"^  rad.;  an  improvement  of  a  factor  of  about  4. 


Figures  43  and  44  display  the  similar  quantity  for  the  other  axis,  i.e., 

PSD  of  LOSy  with  a  disturbance  in  the  x-direction,  both  open  and  closed  loop. 
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The  results  are  similar  also,  with  an  open  loop  LOSy  RMS  error  of  2.0  X  10 
rad.  and  a  closed  loop  value  of  .46  X  10'^  rad.;  again,  a  factor  of  about  4 
improvement. 

Deterministic  Disturbance  Suppression.  The  LOS  responses  to  the  deter¬ 
ministic  disturbances  were  computed  using  a  simulation  of  the  LSS/GTV  facil¬ 
ity  test  scenario.  For  this  reason  the  disturbance  inputs  are  applied  at 
some  finite  "ready"  time  greater  than  zero.  This  simulation  models  the 
structure  continuously  using  a  forth  order  runge-kutta  integration  scheme. 
It  models  the  digital  control  law  discretely  at  the  sample  times  using 
recurrence  algorithms  derived  from  the  z-transforms  of  the  compensators.  The 
structural  model  used  to  produce  the  following  results  included  12  modes  and 
gravity  effects  just  as  the  design  model. 

Fig.  45  shows  LOSx  versus  time  for  a  crew  motion  disturbance  input  in 
the  y-di recti  on  for  both  open  and  closed  loop  cases.  The  crew  motion  distur¬ 
bance  excites  the  bending  of  the  structure  very  little  as  is  evidenced  by  the 
domination  of  LOS^  by  the  pendulum  behavior.  Similarly,  Figure  46  shows  that 
LOSy  due  to  crew  motion  disturbance  in  the  x-direction  is  also  dominated  by 
the  pendulum  mode.  In  both  cases,  the  closed  loop  LOS  curve  shows  the 
increased  damping  and  accompanying  disturbance  rejection  facilitated  by  the 
closed  loop  control.  Figure  47  presents  the  overall  pointing  error  of  Equa¬ 
tion  2-3  for  both  open  and  closed  loop  in  the  face  of  crew  motion  disturbance 
in  both  the  x  and  y  directions  simultaneously.  Again,  the  closed  loop 
response  shows  improved  damping  and  disturbance  rejection. 


The  open  loop  LOS^  error  response  to  the  RCS  thruster  firing  disturbance 
in  the  y-direction  of  Fig.  48  indicates  that  the  RCS  disturbance  excites  the 
bending  of  the  structure  considerably  more  than  the  crew  motion  disturbance. 
The  lack  of  high  frequency  characteristics  in  the  accompanying  closed  loop 
response  indicates  the  excellent  damping  of  the  bending  characteristics. 
This  response  is,  again,  dominated  by  the  pendulum  behavior.  The  LOSy  results 
of  Fig.  49  show  the  damping  of  the  bending  modes  even  more  dramatically  than 
Fig.  48.  Finally,  the  overall  pointing  error  response  of  Fig.  50  due  to  the 
RCS  disturbance  includes  no  surprise  in  its  performance  improvement. 

Conclusions.  The  ASCOT-I  controller  in  conjunction  with  the  ASCOT-I 
structural  model  has  been  shown  to  exhibit  considerable  improvement  in  dis¬ 
turbance  rejection  over  the  open  loop  structural  behavior.  The  improvement 
in  rejection  of  random  disturbances  was  shown  to  be  roughly  a  factor  of  4  in 
each  axis.  The  LOS  responses  to  deterministic  disturbances  indicated  the 
increased  damping  of  both  the  bending  behavior  of  the  structure  and  the 
pendulum  behavior. 

2.1.7  Evaluation  of  which  spacecraft  parameters  affect  the  controller 
performance 

In  order  to  analytically  model  complex  spacecraft,  a  number  of  simpli¬ 
fying  assumptions  must  be  employed  in  the  formulation  of  the  equations  of 
motion.  Typically,  some  of  these  assumptions  involve  the  characterization  of 
the  mass  and  stiffness  properties  of  the  structure,  which  in  turn  directly 
effect  the  spacecraft's  natural  frequencies  and  modeshapes.  Additionally, 
the  characteristics  of  any  actuators  and  transducers  must  be  modeled  in  some 
fashion  such  that  their  dynamics  are  properly  accounted  for  in  the  system 
dynamical  equations. 
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As  far  as  the  RTOP  system  (See  Section  2.2.5)  is  concerned,  the  exact 
orientation  of  the  tip  instrument  package  is  considered  to  be  less  well  known 
than  most  other  parameters.  To  study  the  effects  of  a  change  in  orientation 
of  the  tip  package,  a  model  was  prepared  with  the  tip  package  skewed  off  the 
nominal  position  by  a  30"  angular  displacement  about  an  i  +  j  +  k  vector. 

This  perturbed  model  was  then  used  to  conduct  a  robustness  study  on  the 
rtop  controller  design  described  in  Section  2.2.5.  The  robustness  study  is 
delineated  in  Section  2.1.8. 

2.1.8  Investigation  of  System  Robustness 

The  RTOP  controller  robustness  study  is  very  straightforward.  The  RTOP 
controller  described  in  Section  2.2.5  is  coupled  with  the  perturbed  model  of 
Section  2.1.7  and  the  resulting  closed  loop  system  is  studied. 

The  perturbed  model  does  not  exhibit  significantly  different  in-axis 
behavior  than  that  of  the  RTOP  model.  That  is  to  say,  the  transfer  charac¬ 
teristics  upon  which  the  compensation  was  designed  do  not  change  signifi¬ 
cantly,  and  as  a  reslut  the  system  stability  margins  do  not  change. 

However,  the  cross-axis  characteristics  do  change  significantly  between 
the  RTOP  and  perturbed  models.  Fig.  51  displays  the  open  loop  transfer 
characteristic  from  torque  in  y  (Ty)  to  faceplate  angular  rate  in  x  (©x*) 
the  RTOP  model.  The  same  quantity  for  the  perturbed  model  is  shown  in  Fig. 
52.  Clearly,  the  perturbed  model  exhibits  greater  cross-axis  coupling,  and 
the  robustness  studies  show  that  the  RTOP  controller  is  robust  to  these 
changes.  This  is  largely  due  to  the  fact  that  the  controller  is  not  coupled 
between  axes  and  makes  no  attempt  to  control  one  axis  with  signals  measured 
from  another.  ICAT  allows  the  designer  to  easily  impose  such  conditions  on 
the  control  law  at  the  outset  of  the  design  process. 
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Figure  51 


Cross  Axis  Frequency  Response  from  Torque  in  Y 
to  Faceplate  Angle  in  X  for  the  RTOP  Model . 


2.2  TASK  2  -  TEST  PLANNING  AND  COORDINATION,  TEST  SUPPORT  AND  CONTROL 

Technique  Selection 

2.2.1  Formulation  of  Technical /Management  Program  Development  Plan 

The  Technical /Management  Program  Development  Plan  was  delivered  raid  term 
of  the  contract  period  as  a  stand-alone  document.  It  reviews  the  state  of 
the  ASCOT  effort  in  its  6th  month  and  delineates  the  remainder  of  the  effort 
as  described  in  this  report. 

2.2.2  Modification  of  Ascot- I  to  NASA  Structural  Test  Article 

This  section  provides  a  description  of  the  ASCOT-I  model,  its  origin, 
and  the  motivations  behind  using  it  as  a  test  problem  for  S2D2. 

Introduction.  The  ASCOT-I  model  is  derived  directly  from  the  first  test 
fixture  installed  in  the  LSS  Ground  Test  Verification  (LSS/GTV)  facility  at 
Marshall  Space  Flight  Center  (MSFC).  This  facility  provides  a  test-bed  for 
LSS  type  structures  and  related  control  hardware  and  software  in  both  open 
and  closed  loop  operations.  This  system  and  resulting  structural  model  is 
thought  to  be  a  suitable  test-bed  for  the  S2D2  candidate,  ICAT,  because  of 
the  realism  of  the  system.  Unlike  past  efforts  where  control  strategies  were 
tested  on  LSS  models  derived  from  fictitious  LSS,  the  ASCOT-I  model,  although 
much  simpler,  provides  a  model  of  an  actual  structure  which  has  been  verified 
with  modal  testing.  Also,  reasonable  constraints  on  sensor  and  actuator  per¬ 
formance  and  sample  rate  are  imposed. 

LSS/GTV  Facility  Description.  The  LSS/GTV  facility  and  test  fixture  are 
described  by  the  drawing  of  Fig.  53.  The  first  test  article  will  be  th^’ 
ASTROMAST  beam  as  shown.  The  ASTROMAST  is  extremely  lightweight  (about  5 
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pounds)  and  approximately  45  feet  in  length  and  is  constructed  almost 
entirely  of  S-GLASS. 


The  test  article  is  mounted  on  the  faceplate  of  the  Advanced  Gimbal 
System  (AGS)  engineering  model  which,  along  with  an  additional  torque  actuator 
in  azimuth,  provides  the  control  inputs  for  the  system.  The  azimuth  gimbal 
also  provides  a  means  of  rotating  the  entire  experiment  manually  to  produce 
different  test  scenarios. 

The  gimbal  system  and  associated  electronics  exhibit  similar  dynamic 

behavior  in  each  of  three  axes.  The  transfer  characteristic  from  commanded 

torque  to  applied  torque  is  given  by  a  simple  zero  over  second  order  system 

with  a  simple  underdamped  pole  pair  having  a  natural  frequency  of  100  Hz  and 

a  damping  ratio  of  0.5.  The  transfer  function  is 

Actual  Torque  (2  ir  100)^ 

Gact^  s)  =  - -  =  - - - 1 

Commanded  Torque  s'^  +  (2)(.5)(2  it  100)s  +  (2  n  100)  (2. 2. 2-1) 

The  frequency  response  of  Gact;(s)  is  given  in  Fig.  54. 

The  AGS  is  supported  by  the  Base  Excitation  Table  (BET)  of  the  beam  con¬ 
tainment  structure  which  is  free  to  translate  in  the  horizontal  plane  and 
includes  hydraulic  actuators  to  provide  translational  disturbance  inputs  to 
the  test  fixture.  These  disturbances  can  represent  Astronaut  push-off,  RCS 
(Reaction  Control  System),  thruster  firing,  or  a  broadband  random  disturbance 
such  as  would  be  expected  on  free  flyer  type  spacecraft. 

Six  separately  packaged  inertial  measurement  assemblies  comprise  the 
available  control  system  sensors.  Two  of  the  packages,  containing  three-axis 
translational  accelerometers,  are  identical.  One  is  mounted  on  the  mast  tip, 
and  the  other  on  the  AGS  side  of  the  BET.  Three  other  packages  contain  ATM 
(Apollo  Telescope  Mount)  rate  gyros  and  are  installed  on  the  AGS  faceplate. 
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The  sixth  package,  the  Kearfott  Attitude  Reference  System  (KARS),  is  placed 
at  the  mast  tip  along  with  the  accelerometer  package. 

Each  of  the  sensors  is  assumed  to  have  the  dynamic  characteristic  of  the 

ATM  rate  gyros,  because  this  is  considered  to  be  a  worst  case.  The  transfer 

function  of  a  simple  axis  of  the  ATM  rate  gyros  is 

Measured  Rate  (2  it  40)^ 

Gcen^  s)  =  — — =  —  — 

Actual  Rate  s^  +  (2){.7)(2  ir  40)s  +  (2  n  40)^  (2. 2. 2-2) 

The  frequency  response  of  Gsen(s)  appears  in  Fig.  55.  When  this  transfer 

function  is  used  to  represent  the  accelerometers,  it  represents  Measured 

Acceleration  /  Actual  Acceleration. 

The  signals  from  the  inertial  measurement  instruments  are  read  by  the 
COSMEC-I  data  gathering  and  control  system,  and  processed  50  times  per  second 
according  to  the  control  strategy  under  scrutiny.  The  control  actuator  sig¬ 
nals  are  transmitted  to  the  gimbal  system  as  inputs  to  the  dynamical  system. 

The  COSMEC  1  is  interfaced  to  a  Hewlett  Packard  HP9845C  desktop  computer 
which  stores  data  as  it  is  collected  during  a  test  run,  and  then  provides 
post  experiment  data  reduction  and  display  off  line.  The  controller  inputs 
and  outputs  (measurements  and  commands)  are  recorded  at  each  sample  period  or 
at  some  multiple  of  sample  periods. 

ASCOT-I  Structure  Mathematical  Model.  The  model  of  the  ASCOT -I  struc¬ 
ture  serves  several  purposes  -  each  of  which  are  key  to  the  success  of  the 
analysis.  The  modeshapes  and  frequencies,  which  represent  the  dynamics  of 
the  open-loop  system,  are  extracted  from  the  model  via  the  eigensolution  of 
the  system  equations  of  motion.  These  modeshapes  and  frequencies  are  typi¬ 
cally  the  necessary  input  needed  to  generate  a  closed-loop  control  system. 


An  additional  function  of  the  dynamic  model  is  in  the  simulation  of 
system  response  due  to  a  prescribed  forcing  function.  Clearly,  the  accuracy 
of  the  predicted  response  relies  heavily  on  tiie  accuracy  of  the  mathematical 
model  which  represents  the  physical  system. 

Before  the  procedures  used  to  actually  model  the  ASCOT-I  system  are  dis¬ 
cussed,  the  system  components  considered  the  most  critical,  either  struc¬ 
turally  or  inertially,  will  be  identified.  Referring  to  Fig.  53,  it  is  seen 
that  the  ASCOT-I  system  is  composed  of  six  major  components.  The  six  compo¬ 
nents  are  identified  as  -  1.  the  Base  Excitation  Table  (BET),  2.  roll  gimbal , 
3.  Advanced  Gimbal  System  (AGS)  gimbals,  4.  AGS  Faceplate,  5.  ASTROMAST  beam, 
and  6.  the  tip  instrument  package.  Excluding  the  mass  of  the  BET,  the  mass 
of  the  ASCOT-I  system  is  approximately  310  kg.  A  discussion  of  how  each  of 
the  components  is  modeled  is  presented  in  the  following  paragraphs. 

Modeling  Techniques.  The  main  component  of  the  ASCOT  system  is  the 
ASTROMAST  beam.  For  the  purposes  of  describing  the  dominant  dynamic  charac¬ 
teristics  of  the  system,  the  ASTROMAST  beam  alone  is  considered  flexible. 
All  other  components  of  the  system  are  modeled  as  rigid  elements. 

Over  the  last  several  years,  several  papers  have  concerned  themselves 
with  the  analysis  of  lattice  type  beam  structures  [6,  7].  Although  several 
promising  analysis  techniques  have  been  developed,  it  appears  that  measured 
modal  data  is  presently  more  accurate  than  the  modal  data  obtained  from  the 
analytical  model.  Based  on  this  premise,  a  finite  element  model  of  the 
ASTROMAST  beam  was  assembled  and  "tuned"  to  match  the  measured  modal  fre¬ 
quencies  obtained  from  a  modal  test  performed  at  MSFC.  A  comparison  between 
the  natural  frequencies  presented  by  the  finite  element  model  and  the  mea¬ 
sured  frequencies  is  presented  in  TABLE  2.  As  can  be  seen  in  the  table,  the 


finite  element  model  has  been  tuned  quite  closely  to  the  measured  natural 
frequencies.  As  expected,  the  difference  between  the  predicted  and  measured 
frequencies  tends  to  increase  through  the  higher  frequencies.  This  is  because 
as  the  beam  assumes  more  complicated  displacement  patterns  the  local  non- 
linearities,  which  are  not  explicitly  modeled,  become  more  significant  to  the 
beam  dynamics. 

Since  the  ASCOT-I  system  must  perform  in  a  1-g  environment,  the  effect 
of  Earth's  gravitational  field  on  the  system  dynamics  must  be  evaluated.  The 
gravity  load  on  the  ASCOT-I  system  manifests  itself  mainly  by  altering  the 
frequency  associated  with  the  rigid  body  modes  defined  by  the  hinge  action  of 
the  AGS  gimbals.  These  particular  modes,  shown  in  Fig.  56,  are  essentially 
rigid  pendulum  modes,  with  a  natural  frequency  of  0.15  Hz.  The  value  of  this 
pendulum  frequency  was  obtained  through  the  use  of  Lagrange's  equations 
applied  to  a  1-g  rigid-body  representation  of  the  ASCOT  system. 

The  finite  element  model  of  the  ASCOT  system  was  used  to  produce  the 
system  modes  and  frequencies  (except  for  the  pendulum  frequencies).  The 
natural  frequencies  of  the  system  are  listed  in  TABLE  3,  and  some  representa¬ 
tive  modeshapes  are  plotted  in  Fig.  56.  Except  for  the  pendulum  modes,  the 
"beam"  portion  of  the  ASCOT  system  deflects  in  each  of  the  system  modes  much 
like  a  beam  with  pinned-pinned  boundary  conditions.  This  is  clearly  the 
result  of  the  tip  instrument  mass,  which  tends  to  constrain  the  tip  transla¬ 
tions  since  it  has  a  mass  of  more  than  10  times  the  mass  of  the  beam. 

A  model  of  a  more  complex  structure  which  is  currently  installed  in  the 
LSS/GTV  facility  and  was  obtained  from  MSFC  is  delineated  in  Section  2.2.5. 
This  model  provides  the  basis  for  a  second  and  more  challenging  application 
of  the  S2D2  candidate,  ICAT. 
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2.2.3  Derivation  of  Criteria  for  Evaluating  Candidate  Control  Systems 

Introduction 

To  satisfactorily  ascertain  control  system  performance  in  both  a  com¬ 
parative  and  an  absolute  sense,  a  performance  evaluation  criterion  is  needed 
which, 

1)  uses  the  Line  of  Sight  (LOS)  concept,  thus  preserving  the  spirit  of 
previous  DARPA  sponsored  efforts  in  vibration  suppression  and  point- 
i ng  control , 

2)  is  an  indication  of  vibration  suppression  as  well  as  pointing 
accuracy, 

3)  can  be  measured  and  evaluated  during  structural  tests,  and 

4)  can  be  applied  equally  as  well  to  the  ASCOT  I  and  the  RTOP  (Research 
Technology  Operating  Plan)  models  as  described  in  Sections  2.2.2  and 
2.2.5  respectively. 

The  following  paragraphs  delineate  a  performance  measure  which  meets  the 
goals  set  forth  above.  This  LOS  criterion  is  used  in  the  system  performance 
evaluations  of  Sections  2.1.6  and  2.2.5. 

Evaluation  Criteria.  The  evaluation  criteria  measurements  in  the  y-z 
plane  for  the  ASCOT  I  structure  are  presented  in  Fig.  57.  A  light  source  is 
attached  to  the  tip  instrument  mounting  bracket  such  that  it  points  in  the  z 
direction  and  strikes  a  light  sensitive  position  measurement  system.  With 
reference  to  Fig.  1,  the  point  at  which  the  light  beam  strikes  the  sensor 
with  the  structure  at  rest  (pre-test)  is  called  "A".  Deviation  (denoted  ay 

in  the  y  direction  and  ax  in  the  x  direction)  of  the  light  beam  from  point  A 
would  serve  as  an  acceptable  performance  evaluation  criterion  except  that  it 
would  include  as  error  the  pure  (rigid  body)  translation  of  the  system  induced 
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by  the  x-y  translation  table.  This  is  not  desired  because  ASCOT  is  intended 
to  study  errors  due  to  pointing  and  structural  vibrations.  To  remove  the 
components  of  ax  and  ay  due  to  rigid  body  translations,  the  deviations 
(denoted  tx  and  ty)  of  the  x-y  table  from  its  at  rest  (pre-test)  position  are 
subtracted  from  ax  and  ay  to  yield  measures  of  the  pointing  accuracy  and 
vibration  suppression  performance  of  the  system. 

Using  a  small  angle  approximation,  the  pointing  error  quantities  may  be 
expressed  as 


(2. 2. 3-1) 


for  pointing  error  about  the  x-axis  and 

3x  ■  ^x 


(2. 2. 3-2) 


for  pointing  error  about  the  y-axis.  An  overall  pointing  error  quantity  is 
then  defined  as 


0?^  =  (9^x  +  0y^)l/2 


(2. 2. 3-3) 


Any  or  all  of  the  above  quantities  (Equations  (2. 2. 3-1),  (2. 2. 3-2),  (2. 2. 3-3)) 
may  be  examined  in  determining  the  system  performance.  These  quantities  may 
be  studied  in  terms  of  their  PSD's  and  RMS  values  when  broadband  random 
disturbances  are  applied  to  the  system  or,  in  the  case  of  deterministic 
disturbances,  their  time  functions  may  be  studied. 

Conclusion.  The  evaluation  criteria  delineated  above  satisfy  the  four 
goals  set  forth  in  the  introduction.  In  particular,  they  use  a  line  of  sight 
concept  which  is  an  indication  of  vibration  suppression  and  pointing  accuracy 
but  does  not  involve  pure  translation  of  the  structure.  The  criteria  can 


easily  be  measured  and  evaluated  during  structural  tests  using  relatively 
simple  hardware  which  may  be  affixed  to  the  structure  without  significantly 
changing  the  structural  behavior.  Finally,  these  criteria  may  be  applied  to 
the  RTOP  structure  without  modification,  because  the  same  essential  struc¬ 
tural  elements  compose  both  it  and  the  ASCOT-I  structure. 

2.2.4  Selection  of  Candidate  Control  Systems 

This  section  discusses  briefly  three  well  developed  LSS  control  system 
design  techniques  which  are  considered  candidates  for  application  to  struc¬ 
tures  such  as  the  ASCOT-I  and  RTOP  structural  models. 

Low  Authority  Control/High  Authority  Control  (LAC/HAC).  Lockheed  Mi s- 
sile  and  Space  Corporation's  LAC/HAC  approach  to  LSS  control  system  design  is 
characterized  by  a  two  pronged  attack  of  the  problem.  Low  Authority  Control 
(LAC)  is  used  to  moderately  modify  the  damping  of  the  high  frequency  (i.e., 
less  certain)  modes.  LAC  appears  to  be  readily  implementable  but  the  per¬ 
formance  attributable  to  it  seems  limited. 

High  Authority  Control  (HAC)  is  coupled  with  LAC  to  provide  heavy  damping 
of  the  low  frequency  (typically  better  known)  modes.  HAC  is  implemented 
using  a  Linear  Quadratic  Gaussian  (LQG)  design  which  requires  all  states  of 
the  system  to  be  measured  or,  at  least  estimated.  This  unfortunately  leads 
to  very  high  order  digital  controllers. 

Positivity.  Investigation  of  TRW's  Positivity  technique  has  shown  that 
the  design  is  a  multi step  process  beginning  with  a  predesign  which  modifies 
the  plant  characteristics  so  that  it  is  positive  real.  Then  performance 
specifications  are  satisfied  by  designing  a  feedback  controller  that  is 
positive  real.  This  procedure  assures  that  the  closed  loop  system  will  be 


stable.  However,  because  the  design  is  based  upon  a  truncated  system  model, 
the  actual  system  (proof  model)  Is  not  assured  of  stability. 


Modal  Error  Sensitivity  Suppression  (MESS).  General  Dynamics'  MESS  is  a 
third  candidate  for  LSS  control  system  designs.  The  MESS  design  procedure 
requires  the  separation  of  the  modeled  modes  into  two  groups:  (1)  low  fre¬ 
quency  modes  to  be  controlled  and  (2)  High  frequency  modes  to  be  suppressed. 
An  LQG  problem  is  then  solved  to  produce  the  control  law. 

2.2.5  Comparison  of  Candidate  Control  Designs 

This  report  section  covers  the  effort  performed  under  Task  4. 1.2. 5  as 
delineated  in  Section  IV  of  the  ASCOT  Technical /Management  Program  Develop¬ 
ment  PI  an. 

Introduction.  Application  of  the  S2D2  candidate,  ICAT,  to  the  RTOP 
(Research  Technology  Operating  Plan)  model  of  the  LSS/GTV  facility  fixture  is 
covered  in  this  section.  Discussion  begins  with  a  description  of  the  struc¬ 
ture  and  resulting  structural  model.  Then  the  design  model  is  produced,  and 
the  controller  design  description  ensues.  Once  designed,  the  controller  is 
coupled  to  plant  model  to  evaluate  the  overall  system  performance  with 
respect  to  bandwidth  and  disturbance  rejection. 

RTOP  Structural  Model .  The  RTOP  structure,  so  named  because  it  is  being 
used  in  NASA/MSFC  testing  sponsored  under  a  Research  Technology  Operating 
Plan  (RTOP),  is  an  augmentation  of  the  ASCOT-I  structure.  As  can  be  inferred 
from  Figure  58  ,  the  RTOP  structure  and  resulting  model  contain  all  of  the 
principle  components  of  the  ASCOT-I  structure  plus  an  additional  aluminum 
cruciform  connected  to  the  tip  instrument  package.  The  effect  of  the  cruci¬ 
form  is  to  add  several  lightly  damped  flexible  modes  to  the  structure.  An 
indication  of  the  system  complexity  is  the  spacing  and  number  of  the  system 


natural  frequencies,  which  are  shown  in  TABLE  4.  As  in  the  ASCOT-I  model, 
the  effects  of  gravity  upon  the  structure  are  included;  thus,  the  pendulum 
behavior  appears  rather  than  typical  zero  frequency  rigid  body  modes. 

The  RTOP  model  differs  from  the  ASCOT-I  model  in  that  it  exhibits  light 
cross-axis  coupling,  thus  posing  a  more  apt  application  for  S2D2. 

RTOP  Design  Model.  The  RTOP  design  model  along  with  the  control  system 
topology  is  shown  in  Figure  59.  The  blocks  labeled  I(z)  represent  the  rec¬ 
tangular  integration  carried  out  in  the  strapdown  algorithm  and  are  modeled 
with  the  z-transform,  I(z)  =  Tz/(z-l).  Ghofs)  represents  a  zero  order  hold 
device  and  is  modeled  by  the  transfer  function,  Gho(s)  =  (l-e"Ts)/s. 

The  design  model  incorporates  the  first  40  modes  of  the  RTOP  structural 
model.  This  set  includes  five  modes  with  frequencies  greater  than  the  half 
sample  rate  of  157  rad/sec  (50  Hz)  and  a  largest  model  frequency  of  390 
rad/sec  (62  Hz).  The  modes  above  the  half  sample  rate  are  included  so  that 
their  effects  on  the  sampled  data  system  stability  analysis  can  be  realized. 
The  actuator  and  sensor  models  of  Equations  (2. 2. 2-1)  and  (2. 2. 2-2)  respec¬ 
tively  are  included  as  shown  in  Fig.  59,  and  LOSx  and  LOSy  are  computed  as 
described  in  Section  2.3.3. 

The  measurements  used  in  the  RTOP  design  include  the  three  faceplate 
rotational  rates  (0x,  Qy,  0z)  used  in  the  ASCOT-I  design  and  in  addition,  the 
ASTROMAST  tip  translational  accelerations  in  x  and  y  (x,  y).  The  faceplate 
measurements  are  used  in  the  same  manner  as  in  the  ASCOT-I  design  while  the 
tip  measurements  are  used  to  effect  damping  of  the  pendulum  behavior.  This 
is  discussed  in  greater  detail  in  the  following  section.  The  effector  inputs 
are  the  same  as  in  the  ASCOT-I  model:  gimbal  torque  in  each  of  the  three  axes 


TABLE  4.  NATURAL  FREQUENCIES  OF  RTQP  SYSTEM 


Mode  #  Frequency  (Hz)  Description 


1 

0.00 

2 

0.00 

3 

0.00 

4 

0.14 

5 

0.15 

6 

0.38 

7 

1.10 

8 

1.15 

9 

1.16 

10 

1.21 

11 

1.22 

12 

1.27 

13 

1.28 

14 

1.32 

15 

1.45 

16 

3.02 

17 

3.45 

18 

3.90 

19 

6.55 

20 

6.90 

Rigid  Body  (x-trans) 
Rigid  Body  (y-trans) 
Rigid  Body  (torsion) 
Pendulum  (y-z) 
Pendulum  (x-z) 

1st  torsion 
Tip  reformation 
Tip  reformation 
Tip  reformation 
Tip  reformation 
Tip  reformation 
Tip  reformation 
Tip  reformation 
2nd  BNDG  (x*z) 

2nd  BNDG  (y-z) 

2nd  torsion 
5rd  BNDG  (x-z) 

3rd  BNDG  (y-z) 

4th  BNDG  (x-z) 

4th  BNDG  (y-z) 


Ill 


STRUCTURE 


Figure  59.  RTOP  Design  Model  and  Control  System  Block  Diagram. 


RTOP  Controller  Design.  Examination  of  Figure  59  reveals  that  the 
chosen  RTOP  controller  topology  is  decoupled  between  axes  despite  the  coupled 

nature  of  the  RTOP  structural  model.  (Note  that  the  tip  translation 
in  y  represents  gimbal  rotation  about  x).  The  S2D2  candidate,  ICAT,  can 
handle  not  only  the  cross  axis  coupling  of  the  structure  but  also  the  design 
of  cross-axis  coupled  control  systems.  However,  examination  of  the  cross¬ 
axis  transfer  characteristics  of  the  RTOP  structure  (faceplate  rotational 
rate  about  y  due  to  torque  in  2  for  instance)  reveals  that  closure  of  such 
paths  in  the  controller  would  probably  cause  stability  problems  and/or  fail 
to  improve  system  performance. 

The  x-axis  control  law  architecture  is  given  in  the  block  diagram  of 
Fig.  60.  This  control  strategy  uses  faceplate  angular  rate  (©x*)  ^  ^^mp  the 
flexible  modes,  ASTROMAST  tip  rate  {p’)  to  damp  the  pendulum  mode,  and  face¬ 
plate  position  (^x*)  to  facilitate  pointing.  The  blocks  Dy(z),  D0j^(z),  and 
Eqj^(z)  are  unspecified  2- transforms,  and  Tx  is  the  torque  command  output 
which  is  applied  to  the  gimbal  torquer  characteristic  as  shown  in  Fig.  59. 
The  y-axis  control  law  is  identical  in  form  to  that  of  the  x-axis  as  seen 
from  Fig.  61.  The  z-axis  control  law  is  similar  but  lacks  the  translational 
rate  measurement  as  seen  in  Fig.  62. 

Including  all  three  axes,  the  RTOP  controller  involves  eight  feedback 
paths  and,  therefore,  eight  loops.  Five  of  these,  referred  to  as  "rate 
loops",  are  used  to  augment  the  damping  of  the  flexible  and  pendulum  behavior 
of  the  structure  while  the  remaining  three,  "position  loops"  are  used  to 
facilitate  position  control.  Because  the  loops  are  to  be  closed  one  at  a 
time  but  not  independently,  some  choice  of  closure  order  must  be  made.  Fol¬ 
lowing  the  lead  of  experience,  the  widest  bandwidth  loops  in  each  axis  will 


X-AXIS  CONTROL  LAW 


Figure  62. 


be  closed  first.  The  question  of  which  axis  to  begin  closures  in  is  not 
important  in  the  RTOP  design  because  of  the  light  coupling;  however,  in  a 
system  exhibiting  heavy  coupling,  preference  might  be  given  to  a  particular 
axis. 

TABLE  5  provides  a  summary  of  the  eight  loop  closures  required  along 
with  the  compensation  and  resulting  stability  margins  of  each  loop.  The 
loops  are  closed  in  order  of  the  loop  number,  and  each  time  a  loop  is  closed, 
it  remains  closed  from  that  time  on  in  the  analysis.  The  following  discus¬ 
sion  examines  the  closure  of  loops  1,  4,  and  6;  thus  following  the  progression 
of  the  x-axis  controller.  This  is  not  to  say  that  the  x-axis  controller  is 
designed  Independently.  Quite  to  the  contrary,  examination  of  TABLE  5  shows 
that  the  design  is  being  executed  in  the  three  axes  somewhat  simultaneously. 

The  loop  frequency  response  of  the  compensated  x-axis  faceplate  rate 
loop  (loop  number  one)  appears  in  Fig.  63.  Since  the  compensation  is  simply 
a  gain  of  31dB,  the  uncompensated  response  is  the  same  as  that  of  Fig. 63 
except  shifted  down  in  magnitude  by  31dB.  As  in  the  ASCOT-I  design,  it  is 
desired  that  the  modal  peaks  of  the  compensated  loop  response  extend  beyond 
OdB  in  order  to  produce  damping  of  the  modes,  and  this  is  accomplished  for 
several  of  the  flexible  body  modes.  However,  no  damping  is  added  to  the  pen¬ 
dulum  mode.  This  is  primarily  because  of  the  lag  due  to  the  sensor  and 
actuator  characteristics  (60  degrees  at  the  half  sample  rate)  which  were  not 
modeled  in  the  ASCOT-I  design  where  damping  of  the  pendulum  mode  was  accom¬ 
plished.  Dynamic  compensation  is  not  helpful  in  allowing  greater  gain  in 
this  loop  for  two  reasons:  (1)  lag  due  to  the  sensor  and  actuator  dynamics, 
and  (2)  uncertainty  of  the  modal  frequencies  of  the  higher  frequency  modes. 
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The  polar  plot  of  Fig.  64  shows  that  the  closed  loop  system  with  loop  1 
closed  is  stable  with  stability  margins  of  lOdB  and  55  degrees. 

Closure  of  loops  two  and  three  (faceplate  rate  loops  in  y  and  z)  con¬ 
tinues  in  similar  fashion  to  that  described  above  with  the  cumulative  effects 
of  the  previous  loop  closures  included  in  the  design  of  each  compensator. 
This  brings  the  design  to  loop  four  and,  once  again,  the  x-axis.  The  next 
widest  bandwidth  loop  in  the  x-axis  is  the  tip  translational  rate  loop.  This 
loop  is  included  for  the  express  purpose  of  damping  the  pendulum  behavior. 
One  of  the  obvious  problems  with  damping  the  pendulum  mode  using  the  face¬ 
plate  rate  measurement  is  that  it  is  not  sensed  nearly  as  much  (and  does  not 
peak  nearly  as  high)  as  the  flexible  modes  when  sensed  at  the  faceplate. 
However,  for  small  angles  the  pendulum  mode  is  revealed  at  the  ASTROMAST  tip 
as  pure  translation  in  the  off-axis.  With  the  bending  inodes  damped  by  the 
closures  of  loops  one  through  three,  the  pendulum  mode  peaks  nicely  as  shown 
in  the  compensated  x-axis  tip  translational  rate  loop  frequency  response  of 
Fig.  65.  The  compensation  is  a  gain  of  37 dB  and  negative  in  sign.  The  neg¬ 
ative  polarity  is  necessary  because  positive  translation  in  y  represents 
negative  rotation  about  x  and,  therefore,  the  negative  in  the  compensator 
must  be  included  to  provide  the  correct  (odd)  number  of  negative  signs  around 
the  loop.  The  polar  plot  of  Fig.  66  shows  the  stability  of  the  closed  loop 
system  with  a  gain  margin  of  lOdB  and  a  phase  margin  of  90  degrees. 

Closure  of  the  y-axis  tip  translational  rate  loop  is  similar  to  the  x- 
axis  case  except  that  the  negative  in  the  compensation  is  not  needed  because 
no  negative  occurs  due  to  the  coordinate  system. 

The  final  loop  to  be  closed  in  the  x-axis  is  loop  6,  the  x-axis  position 
loop.  The  design  begins  with  consideration  of  the  open  loop  frequency 


response  from  to  ^x*  first  five  loops  closed.  {See  Fig.  67.)  In 

this  response  both  the  pendulum  mode  and  the  dominant  bending  modes  exhibit 
considerable  damping  because  of  the  closure  of  the  first  five  loops.  As  in 
the  ASCOT-I  design,  the  magnitude  frequency  response  does  not  come  in  at 
-20dB/decade  because  of  the  absense  of  the  rigid  body  mode.  As  before,  an 
integrator  is  specified  as  part  of  the  compensation,  and  the  response  of  Fig. 
68  results  after  inclusion  of  a  53dB  gain  factor.  The  polar  plot  of  Fig.  69 
shows  the  stability  margins  of  6  dB  and  75  degrees. 

The  design  is  completed  with  closure  of  the  y  and  z-axis  position  loops 
in  a  similar  manner  to  that  described  above.  The  resulting  compensation  and 
stability  margins  for  these  and  all  of  the  loops  are  given  in  TABLE  5. 

RTOP  System  Performance.  The  RTOP  system  performance  is  evaluated  in 
the  same  manner  as  the  ASCOT-I  system.  (See  Section  2.1.6.)  All  inputs, 
outputs,  and  criteria  are  the  same,  the  difference  between  the  two  analyses 
being  the  structural  model  and  controller. 

The  first  point  of  consideration  is  the  position  loop  bandwidth.  Figs. 
70,  71  and  72  display  the  closed  position  loop  frequency  responses  for  the  x, 
y  and  z  axes  respectively.  The  system  exhibits  a  closed  position  loop  band¬ 
width  of  .9  rad/sec  in  both  the  x  and  y  axes.  The  lack  of  improvement  in 
bandwidth  in  the  RTOP  system  over  the  ASCOT-I  system,  despite  the  more  com¬ 
plicated  eight  loop  control  system,  may  come  as  a  surprise  to  the  reader.  The 
case  is  that  the  improvement  possible  because  of  the  additional  control  loops 
was  offset  by  the  inclusion  of  sensor  and  actuator  dynamics  which  weren't 
present  in  the  ASCOT-I  system.  The  z-axis  exhibits  a  bandwidth  of  2  rad/ sec 
and  is  considerably  more  behaved  than  the  other  axes  because  of  the  lack  of 
multitudinous  flexible  body  modes. 


Figure  70.  Closed  X-Axis  Position  Loop  Frequency  Response  (RTOP 


In  rejection  of  random  disturbances,  the  RTOP  system,  again,  exhibits 
only  slight  improvement  over  that  of  the  ASCOT-I  system.  Fig.  73  presents 
the  PSD  of  the  LOSx  response  to  a  broadband  random  disturbance  in  the  y-direc- 
tion  with  no  control  system  (open  loop),  and  Fig.  74  displays  the  same  quan¬ 
tity  with  the  control  system  included.  Clearly,  much  of  the  peaking  is 
reduced  in  the  closed  loop  case,  resulting  in  a  reduction  in  RMS  error  from 
1.6  X  10"^  rad.  in  the  open  loop  case  to  3.7  X  10“^  rad.  in  the  closed  loop 
case.  Similar  quantities  appear  for  LOSy  evaluation  in  Figs.  75  and  76  which 
represent  a  reduction  of  RMS  error  from  2.0  X  10"^  rad.  to  3.6  X  10"^  rad. 

Disturbance  rejection  in  response  to  deterministic  disturbances  is  com¬ 
puted  the  same  for  the  RTOP  system  as  for  the  ASCOT-I  system.  The  crew 
motion  disturbance  responses  are  considered  first  by  examination  of  Figs.  77, 
78,  and  79  which  display  the  LOSx,  LOSy,  and  overall  LOS  responses,  respec¬ 
tively,  to  the  crew  motion  disturbance.  In  all  cases  the  crew  motion  distur¬ 
bance  excites  the  bending  modes  very  little.  As  in  the  ASCOT-I  studies,  the 
crew  motion  responses  are  dominated  by  the  pendulum  behavior.  However,  the 
RTOP  model  closed  loop  responses  exhibit  considerably  more  damping  than  their 
counterparts  in  the  ASCOT-I  studies.  (See  Figs.  45,  46,  and  47.) 

The  RCS  thrusters  firing  disturbance  responses  provide  a  better  evalua¬ 
tion  of  the  effects  of  the  RTOP  controller  on  the  flexible  behavior  of  the 
structure.  The  open  loop  LOSx  response  of  Fig.  80  shows  that  the  bending 
modes  of  the  structure  are  excited  by  the  RCS  disturbance.  The  corresponding 
closed  loop  response  shows  that  the  bending  behavior  is  quickly  damped  out; 
however,  both  responses  are  still  dominated  by  the  pendulum  action.  Examina¬ 
tion  of  Figs  81  and  82  reveal  similar  behavior  in  the  LOSy  and  overall  LOS 
responses. 
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Figure  74.  PSD  of  X-Axis  LOS  Closed  Loop  Response  to 
a  Random  Disturbance  in  Y  (RTOP). 
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RTOP  System  Conclusions.  The  S2D2  candidate,  ICAT,  has  been  demonstrated 
through  application  to  a  coupled  3-axis  system  having  40  modes,  3  inputs,  5 
measurements,  and  an  8  path  control  law.  The  system  model  also  includes 
realistic  sensor  and  actuator  dynamics  and  the  effects  of  sampling  upon  sys¬ 
tem  stability. 

The  results  are  promising.  The  closed  loop  system  was  shown  to  be 
stable  and  provide  some  amount  of  disturbance  rejection.  Greater  bandwidth 
and,  hence,  disturbance  rejection  than  in  the  ASCOT-I  system  studies  was  not 
achievable,  despite  the  greater  complexity  of  the  control  system,  because  of 
the  deleterious  effects  of  the  sensor  and  actuator  dynamics.  The  system  did, 
however,  produce  greater  damping  of  the  pendulum  behavior  than  the  ASCOT-I 
controller. 
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3.  CONCLUSIONS 


The  Simplified  Systematic  Digital  Design  {S2D2)  candidate,  ICAT  (One 
Controller  At  a  Time),  has  been  shown  to  be  a  viable  candidate  for  design  of 
control  systems  for  LSS.  The  design  process  has  been  demonstrated  on  a 
decoupled  three-axis  LSS  pointing  system  model  (ASCOT-I)  and  on  a  more 

complicated  coupled  three-axis  system  (RTOP)  of  a  similar  design  having  40 
modes  in  the  design  model.  Reasonable  sensor,  actuator,  and  sample  rate 

constraints  were  included,  and  the  technique  is  equally  useable  with  these 
constraints.  ICAT  also  meets  the  intrinsic  quality  requirements  for  S2D2 
delineated  in  Section  2.1.1.  In  particular,  the  analysis  and  design  are  con¬ 
ducted  in  the  digital  domain  and  relative  stability  of  the  closed  loop  system 

can  be  ensured. 

With  respect  to  the  non-uniform  sampling  rate  prefilter  studies,  it  is 
clear  that  more  theoretical  background  must  be  developed  in  this  area  before 
implementable  filtering  techniques  are  possible.  The  fast  uniform  rate  pre¬ 
filter  concept  is,  on  the  other  hand,  straightforward  and  presently  implemen¬ 
table.  The  only  outstanding  drawback  being  the  possibly  limited  number  of 
applications  where  such  a  prefilter  would  be  both  functional  and  cost  effec¬ 
tive. 

Two  areas  of  suggested  future  effort  stand  out.  First  is  the  use  of  an 
automated,  multi-variable,  frequency  domain  tuning  algorithm  to  optimize  ICAT 
controller  designs  to  certain  specifications,  such  as  bandwidth  and  stability 
margins,  in  an  iterative  manner.  Such  an  algorithm  would  allow  more  precise 
"tuning"  of  stable  designs  produced  using  ICAT. 

Additional  effort  should  be  placed  in  the  non-uniform  sampling  rate 
filter  area,  also.  This  includes  investigation  of  other  schemes  for  choosing 


the  sampling  rate  at  each  sample  period  and  the  possibility  of  digital  com¬ 
pensation  design  using  the  non-uniform  sampling  rate  so  that  no  prefilter 
would  be  necessary. 
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